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In this paper we continue our study of the thermodynamics of large gauge theories 
on compact spaces. We consider toroidal compactifications of pure SU{N) Yang-Mills 
theories and of maximally supersymmetric Yang-Mills theories dimensionally reduced to 
+ 1 or 1 -|- 1 dimensions, and generalizations of such theories where the adjoint fields 
are massive. We describe the phase structure of these theories as a function of the gauge 
coupling, the geometry of the compact space and the mass parameters. In particular, we 
study the behavior of order parameters associated with the holonomy of the gauge field 
around the cycles of the torus. Our methods combine analytic analysis, numerical Monte 
Carlo simulations, and (in the maximally supersymmetric case) information from the dual 
gravitational theories. 
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1. Introduction 

In this paper we study one and two dimensional Euclidean SU{N) Yang-Mills theo- 
ries compactified on a circle and on a torus, respectively. This is relevant in particular for 
studying these theories at finite temperature. In the 't Hooft large limit, we determine 
the phase diagrams of these theories as a function of masses (of adjoint matter fields), 
coupling constants and compactification parameters. The systems we study exhibit rich 
dynamics; in particular, they undergo sharp large phase transitions upon varying pa- 
rameters. These phase transitions are associated with the spontaneous breakdown of the 
Win symmetry W — e^W, where W is the holonomy along a non-contractible cycle of 
the compactification manifold. 

Our analysis employs different techniques in different regimes of parameter space. 
When all matter fields are very massive compared to the scale set by the gauge coupling 
(more precisely, when rn^~'^ ^ A, where m is the mass scale of the matter fields, d is the 
dimension of space-time and A = Qym-^ Hooft coupling), the gauge theories we 

study are weakly coupledS, and may reliably be analyzed in perturbation theory. Already in 
this analytically tractable regime our systems undergo sharp phase transitions and display 
a rich phase structure. Outside this perturbative regime, the phase diagrams we study can 
be constrained by the requirement that they reproduce well known results in special limits. 
The strong coupling behavior of maximally supersymmetric Yang-Mills theories may be 
analyzed using generalizations of the AdS/CFT correspondence. Finally, it is sometimes 
practical to employ Monte Carlo simulations to supplement the information from our other 
techniques. Employing all these methods, we are able to present a reasonably complete 
picture of the phase diagram of the systems we study. 

The investigations reported in this paper are similar in spirit to the recent study of 
Yang-Mills theory on x 5"^ It may be recalled that the phase transitions discovered 

in may be thought of as a weak coupling continuation of the gravitational Hawking- 
Page transition in AdS^ x S^. At least some of the phase transitions we study in this paper 
also have gravitational analogues. For example, as we have already reported in the letter 
1^ (building on the work of 0,^,0), the phase transition in maximally supersymmetric 
Yang-Mills theory on is holographically dual to a Gregory-Lafiamme black hole/black 
string phase transition in type II string theory. See section 7 for other examples. 

^ Recall that the massless gauge field has no dynamical degrees of freedom. 
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The structure of this paper is as follows. In section 2, we discuss various general 
properties of low dimensional gauge theories on tori, and in particular the order parame- 
ters that distinguish the various phases that appear. In section 3, we review the behavior 
of Yang-Mills theory dimensionally reduced to dimensions; this is relevant to the small 
volume limit of the higher dimensional cases. In section 4, we use analytic and numerical 
techniques to map out the phase diagram of pure p + 1-dimensional gauge theory dimen- 
sionally reduced to 0+1 dimensions, and its generalizations with masses for the adjoint 
scalars. In section 5, we study the thermodynamics of 1+1 dimensional Yang-Mills theory 
with massive or massless adjoint scalars on a spatial circle by studying the partition func- 
tion of Euclidean two dimensional Yang-Mills theory on T^. In the case where the scalars 
are very massive, wc obtain a fairly complete picture of the rather rich phase structure 
through analytic analysis, leading to the phase diagrams in figures 14 and 15. In sections 6 
and 7, we repeat the analysis of sections 4 and 5 for the maximally supersymmetric Yang- 
Mills theories in + 1 and 1 + 1 dimensions, respectively. For these theories, knowledge of 
the dual gravitational theories provides additional information about the strong coupling 
behavior of the gauge theory. We find that the supersymmetric theories show rather dif- 
ferent qualitative behavior from the non-supersymmetric theories; for the two dimensional 
case this behavior is summarized in the phase diagrams of figures 18 and 20. 

2. Low Dimensional Yang-Mills Theory on Tori: Generalities 

Massless vector fields have no propagating degrees of freedom in one and two dimen- 
sions. As a consequence, pure Yang-Mills theory in these dimensions is exactly solvable, 
and exhibits relatively tame dynamics. However, the same theory displays rich dynami- 
cal behavior when coupled to matter fields. In this paper we will study large N SU (N) 
Yang-Mills theories coupled to matter fields in the adjoint representation. We begin in 
this section with an overview of the properties of these theories. 

2.1. Quantum mechanical gauge theories (d = 1) 

Different SU{N) quantum mechanical gauge theories with adjoint scalar fields can 
behave in at least three qualitatively distinct ways in the IR. In some theories an attractive 
effective potential (at long range) between the eigenvalues of the scalar fields ensures that 
the vacuum state is normalizable, and that the spectrum is gapped. In other theories, the 
scalar potential has exactly fiat directions and the spectrum is ungapped. In yet other 
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theories the long-range scalar effective potential is repulsive; such theories lack a vacuum 
state. 

Theories with a mass gap (see section 4 below for a detailed study of one class of 
examples) may be thought of as one dimensional analogues of confining theories in d = 4. 
In the 't Hooft large limit with fixed A = Qym-^ 0' spectrum of such theories is 
expected to display a Hagedorn-like growth in the density of states, with a high energy 
density of states p{E) ~ exp{E /Th) for energies which remain finite in the large limit 
(see section 2 of for a simple example). It follows that, upon heating, these theories un- 
dergo a 'deconfinement' phase transition at or below their effective Hagedorn temperature 
Th- In other words, the Euclidean partition function on undergoes a phase transition 
when the circumference of the compactification circle is larger than or equal to the inverse 
Hagedorn temperature. The free energy jumps from 0{1) at low temperatures to 0{N'^) 
at high temperatures (where the theory is weakly coupled), so limTv^oo F{T)/N'^ may be 
viewed as an order parameter for the deconfinement transition. As in four dimensions, 
the expectation value of the Polyakov loop (tr(Pexp(z ^ Aq))) in the Euclidean theory is 
another order parameter for this phase transition i. 

Theories of the second type (those with a well-defined vacuum but an ungapped 
spectrum - see section 6 below for a detailed study of an example) may be thought of 
as the analogues of conformal theories in 4 dimensions. We expect such theories to be 
'deconfined' at all temperatures. If such theories undergo phase transitions as a function 
of temperature (this may or may not happen), the two phases are distinguished by an order 
parameter more sensitive than the expectation value of the Polyakov loop or F{T)/N'^, 
each of which is nonzero in both phases. 

We will not consider theories that lack a vacuum state in this paper. 

2.2. d = 2 gauge theories 

We begin by considering 1-1-1 dimensional large SU{N) Yang-Mills theories, with 
purely adjoint matter content, on IRxtime. Our general expectations for the thermody- 
namics of such systems are analogous to those for quantum mechanical systems reviewed 

^ The Polyakov loop vanishes at large circle radius but is non-zero at small radius (there is a 
subtlety in this statement; see footnote 5 below). Intuitively, the Polyakov loop vanishes at low 
temperatures in a theory with a mass gap because only an infinite number of adjoint photons 
can efi^ectively screen a fundamental charge; the mass gap ensures that such a configuration has 
infinite energy. See section 5.7 of Q for a more careful and detailed discussion. 
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in the previous subsection. Gapped theories are expected to have a Hagedorn growth in 
their density of states for energies which remain finite in the large N Umit (see for 
examples). These theories are expected to undergo a deconfining phase transition at or 
below the Hagedorn temperature. The free energy is expected to be of order 0{N'^) (and 
the Polyakov loop nonzero) at high temperatures, while the free energy is expected to be 
(9(1) (and the Polyakov loop vanishes) at low temperatures. On the other hand the free 
energy in ungapped theories is expected to be 0{N'^) (and the Polyakov loop expected to 
be nonzero) at all temperatures i. 

However, field theories in d = 2 have additional structure absent in c? = 1. It is 
possible, and rather natural, to study the thermodynamics of such theories compactified 
on a spatial circle. In Euclidean space this thermodynamics is described by the partition 
function of the corresponding Yang-Mills theory on a T^. This opens the possibility for a 
much richer phase structure in these models, as we explain in the rest of this subsection. 

Recall that we are interested in SU{N) theories which only have fields in the adjoint 
representation. For such theories, the gauge group is really SU{N) /'Eij^ , since the gauge 
transformations in the 'K n center of the gauge group act trivially. Whenever such a theory 
is compactified on a torus, it possesses a (!ZAr)o x (!^Ar)i global symmetry. This global 
symmetry is generated by gauge transformations G that are not periodic on the torus but 
instead obey G{x'^ +peo+qei) = (3^G[x^), where cq and ei are the fundamental cycles of 
the torus (along the and 1 directions) and a and (3 are both gauge transformations in the 
center 'K ^ . All local gauge-invariant operators are uncharged under this global symmetry, 
but fundamental Wilson lines that wrap around a (p, q) cycle before closing carry charge 
(p, q) . Thus, non-zero expectation values for such Wilson loops break this global symmetry. 
In particular, an expectation value for the Wilson loop § spontaneously breaks the 
(!^7v)^i symmetryil, where = 0, 1. 

^ This is true generically, but it is not true in some supersymmetric examples in which the low 
temperature free energy is of order N. 

^ Where W^{x) = j^{tT:{Pe ^ ^^)) around the circle in the ^ direction, for /i = 0, 1. 

^ Actually we should be more careful here - in fact there is no symmetry breaking at finite 
volume. The closest one can come to a symmetry breaking is to have A'^ different saddle 
points (related by the W^n symmetry) dominating the path integral, but the path integral sums 
over all of them. Thus, there will not really be an expectation value for W^, but just for |VK/ip 
(for example). For the most part we will leave this issue implicit in our discussions. 
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As we will see later in this paper, in the decompactification limit i?i — > oo the sym- 
metry (^jv)i is never broken so Wi always vanishes. The Polyakov loop Wq is the sole 
order parameter for the system in this limit, making contact with the discussion at the 
beginning of this subsection. At finite Ri, Wq and Wi are both nontrivial order parame- 
ters, allowing for intricate two dimensional phase diagrams with the four possible phases 
separated by phase transition lines. Of course, such a phase diagram in a system that is 
ungapped in the decompactification {Ri cx)) limit (see section 7 for an example) has 
qualitative differences from its counterpart in a theory that is gapped in the same limit 
(see section 5 for examples). 

In the large limit, instead of just considering the Wilson loop in the fundamental 
representation, it is often usefnl to study the full holonomy matrix U,{1) = Pe' f ^- . 
which is a unitary matrix whose trace (in some representation) gives the Wilson loop 
(in that representation). The set of eigenvalues of this matrix (which live on the unit 
circle) is gauge- invariant, and in the large limit their distribution along the unit circle 
is some continuous function. The (^at)^ global symmetry described above shifts the 
phases of all the eigenvalues of t/^ by an angle of 2tt/N. In the ^jv-symmetric phase the 
large N eigenvalue distribution is constant, while otherwise it is generally maximized at 
some particular value, spontaneously breaking the % n symmetry. Detailed analysis of the 
eigenvalue distributions sometimes permits a sharp distinction between two phases with 
the same symmetry breaking pattern, allowing for the possibility of even more intricate 
phase diagrams (interpolating between a larger number of phases) than those described 
above. 

2.3. On dimensional reduction 

One of the methods we will use in our analysis is dimensional reduction. Whenever 
one of the circles in the problem is very small compared to the other scales, the theory may 
be approximated by a theory in one dimension less. So, the functional integral of a (d+ 1)- 
dimensional gauge theory will approximate, in the limit of a very small circumference 
Rd+i for one of the compact circles, that of a d-dimensional gauge theory, with a coupling 
constant = \d+\/ Rd+i- When the small circle has periodic boundary conditions for 
fermions, the lower-dimensional gauge theory has precisely the same field content as the 
original theory, except that the zero mode of one of the components of the vector field 
becomes a scalar field. When the fermions have anti-periodic boundary conditions on the 
small circle (as for thermal boundary conditions on a temporal circle) they are projected 
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out, leaving only the bosonic fields. Naively, the lower dimensional gauge theory is valid 
whenever the Kaluza-Klein scale 1/Rd+i is much larger than the other compactification 
scales 1/Ri and T, and than the dynamical scale Kd = A^*^^ the latter condition is the 
same as Rd+i < >^d+i'^^~'^^ = 1/Ad+i. 

This argument is naive since an SU{N) gauge theory on a torus actually has excita- 
tions whose classical energies scale as l/{NRd+i), coming from configurations for which 
gauge-invariant sets of eigenvalues are periodic only up to a permutation (these configu- 
rations play a major role in M(atrix) theory). This scale goes to zero in the 't Hooft large 
limit, so naively we can never ignore all the KK modes in this limit. The importance 
or otherwise of these modes (and so the validity of dimensional reduction) is a dynamical 
issue determined by the saddle point eigenvalue distribution of the Wilson line operator 
around the corresponding circle. The light modes described earlier in this paragraph are 
present only when the Wilson line eigenvalues are spread over the circle (the 'shift' distri- 
bution in M(atrix) theory language). When the eigenvalues of the Wilson line operator are 
sufficiently clumped, these light modes are absent and dimensional reduction is justified. 
In this paper we will only use dimensional reduction when this criterion is met. 



3. Bosonic Matrix Integrals {d = 0) 



In this section we briefiy review the behavior of bosonic SU{N) matrix integrals; we 
will find the results of this section useful in our analysis below. 
Consider the matrix integralil 



V^i exp 



2Af 



-Tr 



(3.1) 



where z, j = 1, . . . ,p, [p > 2) and the $i are N x N Hermitean bosonic matrices, in the 't 
Hooft large limit in which is taken to infinity with constant Aq. Redefining variables 



(pi, we find 



V{pi exp 



-A^Tr 



* J i<j 



^ 2 



(3.2) 



^ Note that in zero dimensions the only difference between a gauged matrix integral and a 
non-gauged integral is the volume of SU{N), so we will not need to distinguish the two. 
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Let us first study the limit Xo/mf ^ (for alH). In tliis limit the integral factorizes into a 
product of p identical integrals, each of which is easily solved by saddle points. The saddle 
point eigenvalue distribution for (say) fi is given by the Wigner semi-circle law, 

2 



p{x) 



1 - X" 



TT 



and we have 



l(try$2) = ^yi 



2 

mf 



(3.3) 



(3.4) 



In the simple limit considered above, suitably normalized gauge invariant expectation 
values (like ;^Tr^^$^) are completely determined by a saddle point with sharp edges, 
and are independent of A^. This follows from rather general considerations (factorization 
and 't Hooft scaling), and so generalizes to most of the systems we study in this paper. 



<lTr 
N 



i 

3r 



2.5 

2 

1.5 

1 

0.5 



6 



10 



Figure 1: Graph showing (;^Tr^j$f) for the massless matrix integral 
( p.l|) with Ao = 1 for various values of p from 3 to 10. The red points 
correspond to = 10, the blue to = 20. 




Figure 2: Plots showing the distribution of eigenvalues for the massless 
matrix integral with Aq = 1 for p = 3,5,10 (left to right), generated by 
Monte Carlo simulations for = 20. 
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We now turn to the opposite, massless limit of (|3.1|) . It follows from a simple rescaling 

(similar to the one we used above) that the eigenvalues are localized on a length scale 

1 

a = KXq at the corresponding saddle point; Monte Carlo simulations demonstrate that 
K is of order unity (see Appendix A.l for more details and references), and the values of 
(-^Tr $1) are shown in figure 1 for various values of p. Sample eigenvalue distributions 
which demonstrate a qualitatively similar form to (^]^), including its sharp edge, are shown 
in figure 2. 

It natural to guess that ( |3.1| ) is dominated by a saddle point in which the eigenvalues 
of the scalars are sharply localized for all values of nii and Aq. The localization length 
scale varies smoothly from approximately \/Xo/ ^ at large m to approximately Aq at small 
m. This expectation is easily verified (at sample values of parameters) by a Monte Carlo 
simulation. 



4. Bosonic Gauge Theories in One Dimension 

In this section we study the quantum mechanical SU{N) Yang-Mills theory coupled 
to p adjoint scalar fields (p > 2), whose Euclidean action is given by 

where -Dq is a covariant derivative involving the non-dynamical gauge field Aq, we assumed 
equal masses for simplicity, and z, j = 1, ■ ■ ■ This theory is believed to possess a unique 
normalizable vacuum state, and to have a mass gap in its spectrum, at all values of M 
including M = 0. We will study the Euclidean partition function of this theory on a 
circle of circumference i2, which is the same as the thermal partition function of the 
quantum mechanical system at temperature T = 1/R. This theory is characterized by two 

dimensionless parameters, which we take to be the inverse radius in units of the coupling, 

~ i i 

t = l/(i?Af ), and the mass in units of the coupling, m = M/Xl. In this section we will 
determine the phase diagram of (|4.1|) as a function of these parameters. 

Two asymptotic regions of this phase diagram are amenable to analytic analysis. First, 
in the limit of large t, ( [4.1|) effectively reduces to the matrix integral ( p.l|) (as discussed in 
section 2); from the analysis of the previous section it then follows that the eigenvalues of 
the hoJomy maw J Uo ^ Pe'f^' a.e clumped at all values of m at large I 

^ By an abuse of notation, we will sometimes refer to these as the eigenvalues of the Polyakov 
loop operator Wo = jjtr{Uo). 
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At large m, (f4.1| ) may be analyzed in perturbation theory. In the next subsection we 
demonstrate that, in this limit, our system undergoes a first order deconfinement transition 
at tc ~ m/ln(p), and we determine the first correction (in 1/m) to this phase transition 
curve. For t > tc, the eigenvalues of the Polyakov loop operator are clumped and the 
Win invariance is broken. For t < tc the eigenvalues of the Polyakov loop operator are 
uniformly distributed and '^n invariance is restored. 

The analysis of the two limits described above (and especially the analytic form of 
the large m phase transition curve given in the next subsection) suggests that the phase 
diagram of ( |4.1| ) takes the form shown in figures 3 and 4, with the phase boundary extending 
to some fixed temperature for nonzero mass. In sections 4.2 and 4.3 below, we will use 
Monte Carlo simulations to plot out this phase diagram in detail, confirming this prediction 
and giving a quantitative picture of the phase diagram for all values of the mass. 

In section 4.4, we discuss the phase diagram for a closely related theory in which one 
scalar remains massless while the rest have mass M; this theory is relevant to the high 
temperature limit of 1+1 dimensional gauge theories with massive adjoint scalars (to be 
discussed in section 5). We find similar qualitative behavior, except that the transition 
temperature for large mass becomes t ~ m/ln(m). 

4.1. The large mass limit 



When the masses in (4J.) are sufficiently large (m ^ 1), (|4.1|) is weakly coupled at 



all temperatures, and the thermodynamic behavior can be studied in perturbation theory. 
As discussed in [^0, it is straightforward to integrate out all massive degrees of freedom 
in the weakly coupled limit to obtain an effective action in terms of the sole light degree 
of freedom, namely the Wilson line of the gauge field about the circle 

U^Pe'€'''^°. (4.2) 

The partition function then takes the form of a unitary matrix model, 

Z = I DUe~^-'ii^^\ (4.3) 



where 

Seff{U)=pY, — tr(t/'^)tr([/t-) +0(l/m3), (4.4) 



n 

n=l 



with X = e~^^ = e~"^/*. In the weak couphng (m oo) hmit, this matrix model 
(for p > 1) undergoes a large phase transition as a function of t, from a low t phase 
dominated by a saddle point in which the eigenvalues of U are distributed uniformly around 
the circle, to a high t phase in which the eigenvalues of U are clumped in the saddle-point 
configuration. 

In the strict m — oo limit, the phase transition occurs at the point t = m/ln(p), 
and it is (weakly) of first order. However, as discussed in section 6 of 0, the nature of 
the phase transition at large but finite m depends on higher order terms in the effective 
potential, arising from two and three-loop diagrams. In appendix C, we compute the 
effective action ([4.4|) to three-loop order, and find that after integrating out all scalars 
as well as all n > 1 Fourier modes of the eigenvalue distribution function (these are all 
massive near the transition), we are left with an effective action for the lowest Fourier 
mode Ui = j^tr{U), of the form 

SeffM = N^{ml{x, l/m^)\ui\^ + ^6(x, l/m^)\uif + 0{l/m^^)), (4.5) 

where the leading order expressions for ml and b are given in appendix C. Since we find 
6 < at the transition temperature where ml switches sign, this effective action describes 
a phase transition which continues to be first order also for large finite m (as at m — oo) 
(see 0), and occurs at the critical temperature 
1 (p-l)(2p+l) 



tc =m/ ln{p) + 



Am^ p ln{p) 

^ ~ ^) (ln(p)(21p2 + Qp) + 48/ + 20p^ - I7p^ - 20p - 4) + 0{l/m^ 



128m^ p^ lii{p) 

(4.6) 

which is slightly below the Hagedorn temperature of this theory. 

In summary, when m ^ p^ (so that perturbation theory is reliable), ( [4.1| ) undergoes 
a first order phase transition at tc given by (|4.6|) . For t < tc the Polyakov line eigenvalue 
distribution is uniform. For t > tc ^ m/ln(p), the eigenvalues of the Wilson line are 
clumped. 

It is tempting to extrapolate these results beyond the validity of perturbation theory. 
As displayed in figure 3, plotted for p = 4, the phase transition line asymptotes to the line 
t = m/ ln{p) for large m, and gradually rises above it as m decreases, until perturbation 
theory is no longer valid (at roughly m ~ p^). A naive extrapolation of our weak coupling 
results suggests that the phase transition curve will hit the vertical axis (m = 0) at a finite 
temperature. 

In the rest of this section we will use Monte Carlo techniques to demonstrate that this 
guess is indeed correct. 
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Figure 3: A plot showing tiie prediction of a first order phase transition 
at large m» p^/^ from perturbation theory. The small circle limit suggests 
localization of the Wilson loop eigenvalues at large t, for fixed m. Then, 
extrapolating the blue phase boundary would naively indicate that it meets 
the t axis. 
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Figure 4: A plot showing the conjectured phase boundary translated into 
the Ai, plane, choosing units with R= 1. The Monte Carlo data presented 
below confirms this picture. 

4-2. Results for m = from Monte Carlo simulations 

The bosonic one-dimensional model is rather simple to implement numerically 
using elementary Monte Carlo methods. In this section we will discuss the numerical 
results for the behavior of the Wilson line eigenvalues in the cases p = 2, 4, 9, characterized 
by Un = -^|tr(C/"^)| (which vanish for all non-zero values of n in the uniform distribution). 
For these diverse values of p we will see that they all have the same qualitative behavior. 
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and we will link the results to the analytic limits discussed above. Details of the method 
may be found in Appendix E. 

In the massless limit, ( |4.1| ) is simply the dimensional reduction of the pure SU{N) 
gauge theory in (p + 1) dimensions to one dimension. In this limit ( [4.1|) has a single 
dimensionless parameter Ai-R^ = As we noted above, at large t this system reduces 

to the matrix integral of the previous section (with masses set to zero), and so it lies in a 
phase with a clumped distribution of eigenvalues for the Polyakov line. We now describe 
results of a Monte Carlo analysis we have performed on this theory. 
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Figure 5: Figure of (ui) as a function of Ai for p = 9 with various values 
of A^. We choose units with R = 1, and thus Ai = We see that as A^ 

increases the points to the right of Ai ~ 1.4 appear to decrease, consistent 
with 1/A^ scaling. To the left of this value the points appear to tend to a 
limiting curve. This indicates a sharp discontinuity in ui at infinite A^, with 
Ui being an order parameter. Statistical error bars are comparable to the 
point sizes. 



In figures 5 and 6 we plot the Monte Carlo results for (ui) and (^2) in the theory 
with p = 9, for various values of A^. These results were previously reported in 0; as we 
have reported there, our data strongly suggests a large A^ transition at ^ 1.4. As is 
apparent from figures 5 and 6, the low t (or large Ai) phase has ui = and U2 = (more 
precisely, ui and U2 are of order 1/A^ for finite A^). Clearly, ui and U2 are nonzero (and 
hence the eigenvalue distribution is non-uniform) in the high temperature phase. All this 
is perfectly in line with the conjecture made in the previous subsection. 

As discussed in [P], we have not been able to clearly resolve the order of the phase 
transition studied in this subsection. Unfortunately, the phase transition of our system 
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12 3 4 

Figure 6: Figure of (^2) as a function of Ai for p — 9 witti various values 
of A^. Again a discontinuous behavior is indicated in the large limit. As 
above, we choose R — 1 units. 

appears to lie very near the boundary between first and second order behavior (see 0). 
In this situation it is difficult to numerically distinguish between the two reasonable possi- 
bilities (see Appendix E), which are either a first order phase transition or a second order 
phase transition followed by another, third order phase transition 0. However, it is clear 
from the data that if the second scenario is correct then the two phase transitions must be 
very close together. 

We have repeated the analysis of this section aX p = 2,4. We find qualitatively 
identical behavior to that reported in the previous paragraphs, but this time with the 
phase transition at XcritR^ = l/^c — O-^, 0.9, respectively. 

4-3. Monte Carlo simulations with non-zero mass 

We have also performed Monte Carlo simulations of ( [4.1| ), as a function of t, at various 
different values of m. The results of these simulations (which we present in this subsection) 
smoothly interpolate between the analytic prediction of subsection 4.1 and the m = 
results of the previous subsection, allowing us to fill out the phase diagram of (^]l|) ^^'^ ^° 
confirm the guess displayed in figure 4. 

In figures 7 and 8 we present three dimensional plots of {ui) and (^2) as functions of 
XiR^ = 1/? and MR = m/t, for the case p = An. An examination of figures 7 and 8 
reveals that the m, t plane may be divided into two regions. In the first (uniform) region ui 

^ These graphs are best viewed in colour. 



U2 

1 



. 8 



. 6 



0.4 



0.2 



13 



p=4, surfaces of u 



1n 




Figure 7: Figure of {ui) in the Ai,M plane (in R — 1 units). The red 
(top) surface is for = 6, blue is = 12 and green (bottom) = 20. We 
see a transition line bounding the origin where outside the line ui decreases 
with increasing (approximately consistent with 1/A^ scaling), and inside 
the line a limit is approached. The transition line appears to connect the 
massless transition to the large m transition. The black mesh of the surfaces 
gives the sampling density. Thus, the (small) statistical error can be seen as 
the 'roughness' of the surface. 



p=4, surfaces of u 




Figure 8: Figure of (^2) in the Ai, M plane, employing the same conventions 
as in the previous figure. 



and U2 decrease as A^ increases; this decrease is approximately fit by a 1/A^ decay to 
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In the second (clumped) region ui and U2 asymptote to fixed, m and t dependent, nonzero 
values. These two regions are divided by a transition line. As in the case of m = 0, ui 
appears to jump disco ntinuously (or at least with a very sharp slope) across this transition 
line in the limit N —>■ oo. 




0.1 0.2 0.3 0.4 0.5 

X 



Figure 9: Contours of A^^'^'^(ui) for the p = 2 theory, in the Ai,M plane 
(in R = 1 units). The green region is predicted to be in the uniform eigenvalue 
phase, and the red region in the localized phase. The blue curve gives the 
large m perturbative prediction, where the perturbative expansion is valid. 
We clearly see that the phase boundary connects the two axes. 

In figures 9, 10, 11 we plot this transition line as a function of m and t (see below for 
more details on how exactly this transition line was obtained from our data) . In each case 
this line interpolates smoothly between the analytic prediction of subsection 4.1 and the 
results of the previous subsection. As in the previous subsection, our numerics are unable 
to definitively establish the order of the transition; however they are certainly consistent 
with the simplest conjecture, which is that the transition (which the computations of 
subsection 4.1 have established to be of first order at large m) remains of first order for all 
values of m. 

To end this subsection we give a precise definition of the transition lines plotted in 
figures 9, 10, 11. We characterize the position of the transition by computing the difference 
between the values of (ui) for different values of A^. 'Outside' the phase transition boundary 
ui decreases with increasing A, while 'inside' it stays roughly constant (with a small 
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X 



Figure 10: Contours of A^^'^'^(iii) for the p — 4 theory, with the same 
conventions as above. 




Figure 11: Contours of A^'^^(w,i) for the p — 9 theory, with the same 
conventions as above. 

increase). We may estimate the phase boundary as the line where these differences are 
zero, namely the function ui is neither increasing nor decreasing with increase of A^. Of 
course, this estimate to the phase transition line becomes more accurate as we increase A^. 
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Figure 12: Figure of A^'^^{ui) (blue) and A^'^''^^{ui) (red) in the Ai,M 
plane for the p = 4 theory (taking R = 1 units). We take the zero contour to 
measure the position of the phase boundary. 

To illustrate our procedure, in figure 12 we plot A^'^*^(wi) = (ui)^^^^^ — (ui)^^^"^^^ 
and A^^'^°('Ui) = {ui}^^^^"^^ — {ui)^^^'^^\ Both surfaces in the plot give a consistent phase 
boundary location. 

4.4- One massless and p massive scalar s 

Before concluding this section, we briefly consider a generalization of ( |4.1| ) for which 
one of the scalar s (which we call (f)) remains at zero mass, while the other p scalars are 
massive. This arises in studying the high temperature limit of a 1+1 dimensional gauge 
theory with massive scalars, for which the temporal component of the gauge field leads to 
a massless scalar in the dimensionally reduced theory. 

The theory behaves in a qualitatively similar manner to the massive theory considered 
above. For m 0, the two theories are the same, while at large m, we will see that both 
theories have a phase boundary at a temperature which goes to infinity as m goes to 
infinity, though with different dependence on m. 

In this limit, it is again appropriate to integrate out the massive scalars at one loop 
to generate an effective potential for cf) and the holonomy matrix U . The (j) independent 
terms will be given by ( [4.4| ) as before, with corrections negligible as long as m ^ p^/^. At 
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very large m, the most important (p-dependent term will be a mass terml 

Seff = I dt^M^). (4.7) 



2M 

The leading (j)'^ interaction comes with a coefficient g"^ ~ \\/NM^^ and can be ignored 
relative to the quadratic terms when integrating out (j) as long as 

m^/(4Ar)^>l ^ m>l. (4.8) 

There will be additional terms involving both cj) and t/, but these are nonlocal and will 
be suppressed exponentially by powers of e~^^ for large R. We should be able to ignore 
these if 

i2M>l =^ t<^m. (4.9) 

When all these conditions are satisfied, we may integrate out the scalar (p at one loop 
(using only the kinetic term and mass term) to obtain additional terms in the effective 
action for U. The result (including terms we have already from integrating out the scalars 
of mass M) is 

S,ff{U) = J2 -b^" + ^"') tr([/")tr(t/^"), (4.10) 

n=l ^ 

where x = e~^^ = e~"^/* as before and e is the ratio of the induced (p mass to the X 
masses, 

e=^«l. (4.11) 
This theory has a deconfinement phase transition when 

px + x"^ = 1 
=^ px + eln(x) = 

=^ X 



^(In(^) + C(ln(ln(^)))) (4-12) 



p e e 

~ 2 m f /Inflnfm)) 
t = \l + 0' ^ ^ 



3 ln(m) \ \ ln(m 



Note that at this temperature, ( |4.8| ) and ( |4.9| ) are both satisfied, so our analysis is self 
consistent. Thus, we conclude that at large m, the phase transition temperature goes as 
t = |m/ ln(m). 



^ Here, cj) is normalized to have a kinetic term J c?ttr(i(Z)o<; 
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5. Two Dimensional Bosonic Gauge Theories on 

In this section we study the two dimensional large SU{N) Yang-Mills theory coupled 
to p Hermitian scalar fields in the adjoint representation, 

S = ^Jd'xtr l^f, + J2 [{d.'^' - ^ [A,. + - Y: [^^ $1 '} , (5.1) 

with /, J = 1, - ■ ■ ,p. Viewed as a 1 + 1 dimensional quantum field theory, this theory 
is believed to possess a unique normalizable vacuum state and a Hagedorn growth in its 
density of states (see for the special case p = 1). We will study the partition function 
of the Euclidean theory on a rectangular torus with radii Ri and R2, as a function of the 
geometry and masses (in units where the coupling is fixed) . We will investigate this system 
in two separate limits; first when all scalars are massless, and second when all scalars have 
large mass. 

In this case, the geometry has two non-contractible cycles, so any saddle point con- 
figuration will be characterized by two holonomy matrices which we call U and V. At 
each point in parameter space, we will ask whether the eigenvalues of these matrices in the 
saddle point configuration are clumped or unclumped (uniformly distributed on the unit 
circle), and so obtain a phase diagram containing four possible phases. The phase diagram 
will turn out to be similar in the massless and large mass limits, and so we believe it is 
likely to have the same structure for any value of the mass. 

5.1. The massless theory 

We first consider the theory ( |5.1|) in the limit where the scalar masses are zero. We 
wish to understand the phase diagram of this theory on a rectangular torus, as a function of 
its two dimensionless parameters, which we can choose to be ri = Ri\/X and r2 = R2y/X, 
the circumferences of the spatial and temporal circles in units where the coupling is set 
to one. Since there is nothing to distinguish the two circles, the phase diagram must be 
symmetric under the exchange ri ^ r2. 

Let us first review what we know about the theory in the limit ri — > 00 keeping r2 fixed; 
this limit gives the decompactified system at a temperature l/i?2- At low temperatures 
this system is believed to confine, with a Hagedorn tower of glueball states (see [|] for a 
proof for the case p = 1). As the temperature is raised, the system undergoes a sharp 
deconfinement transition at a temperature of order V^. In Euclidean space the order 
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parameter for this phase transition is the expectation value of the Wilson line around 
the X2 circle; the eigenvalues of this Wilson line are clumped at high temperatures but 
uniformly smeared over the circle at low temperatures. Thus, for large ri, we expect a 
deconfinement transition line which asymptotes to a constant r2 of order unity. Exchanging 
the two circles, we must then have also a phase transition line for large r2 asymptoting to 
constant ri, across which the eigenvalues of the Wilson line in the xi direction clump as 
we decrease ri. 

We now turn to the opposite limit where the spatial circle i?i is very small. When 

1 /3 

the Kaluza-Klein scale 1/Ri is much larger than the scale set by the effective one- 
dimensional coupling Ai = A/-R1 and also much larger than the temperature I/R2, the 
theory will reduce to the one-dimensional theory considered in the previous section with 
p + 1 massless scalars. The parameter t of the previous section is given in terms of the 
two-dimensional parameters by 

1 

~ 1 r 3 

t= r = — • (5.2) 

i?2Af ^2 

Thus, in the regime ri ^ 1 and ri -C r2 we expect a phase transition along the curve 

r2 = (5.3) 

where t^'^^^ is the critical value of t in the one-dimensional theory with p + 1 massless 
scalars. The eigenvalues of U (the holonomy around R2) will be clumped or unclumped 
for values of r2 below or above this line (again assuming that ri <^ 1 and ri <^ r2). 

Again, we may swap the role of the two circles and conclude that there is an additional 
phase transition along the curve 

r2 = {Vf+^^frl (5.4) 

such that the eigenvalues of V (the holonomy around Ri) are clumped/unclumped 
above/below the curve, in the regime of validity r2 <^ ri and r2 <S 1. 

Note that these results rely on the numerics of the previous section, however, in the 
regime 

rf > r2 > (5.5) 
20 
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Figure 13: Known aspects of the phase diagram of massless bosonic SU{N) 
gauge theories with three possible completions. U is the expectation value 
of the holonomy in the X2 direction, and V is the expectation value of the 
holonomy in the Xi direction. 

all non-zero modes are weakly coupled and can be integrated out at one loop (see appendix 
A) so we are able to verify analytically that both U and V are clumped in this regime.ll^l 

Putting all this together, we conclude that the phase diagram of ( |5.1| ) for m ^ 1 takes 
the form shown in figure 13. We have sketched three possible simple completions of this 
phase diagram, however, our current lack of understanding of the dynamics at ri,r2 ~ 1 
prevents us from distinguishing between these three possibilities. 

The phase transitions in figure 13 are driven by the strongly coupled dynamics of a 
field theories' worth of degrees of freedom, and so are difficult to analytically control. In 
the following subsections, we will see that as for the one-dimensional theories considered in 
section 4, analytic analysis becomes possible in the opposite limit when the scalar masses 
become large, m ^ 1, where we take mi = M and define the dimensionless mass parameter 
m = M/VX. 

Naively, we might expect that the theory should admit a perturbative analysis as long as the 
torus area is small in units of the coupling. In fact, as in finite temperature computations in 3+1 
dimensional theories, this naive expectation is modified by infrared effects (see for a discussion 
of this point). The correct regime of validity may be obtained by requiring that the Ri circle 
should be small in units of the one-dimensional coupling X/R2 and that the R2 circle should be 
small in units of the one-dimensional coupling X/Ri, leading to ( ^.51) . 
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5.2. M = oo 



First, consider the strict limit M — oo. In this hmit all scalar fields may simply be 



set to zero and ( |5.1j ) reduces to pure two dimensional Yang-Mills theory on a torus, an 
exactly solvable system. Over 30 years ago 0, Migdal rewrote the partition function of 
two dimensional Yang-Mills theory on a torus in terms of an integral over the two unitary 
holonomy matrices 



Zyrr,= I DU DV Y^dRC-^^'^'^hRiUVU-W-^), 

R 



where 



A = XR1R2 = rir2. 



(5.6) 



(5.7) 



Here, the sum over R runs over all representations of the SU {N) gauge group, and dn and 
C2{R) are the dimension and the quadratic Casimir of the representation R (see Appendix 
D for a brief derivation). It is not difficult to integrate out one of the two matrices - say, 
V - in (|5.6|) , to obtain an effective action for U , 



■'ym 



duj: 



e 2iv 



C2(R) 



XR 



{U)xr{U^) 



(5.8) 



R 



which may be evaluated to obtain [p!0|JTl] , p!2| , p!3[ 



e 2iv 



^C2{R) 



(5.9) 



R 



For our analysis below, it will be useful to note (see Appendix D) that at large A^, 
) may be written as 



1 



Zym= / I^t/exp Q]]-(-e-^'^ + 2e-^)tr(t/")tr(t/-'^) 



(5.10) 



Equation (|5.10|) may be rewritten as an integral over the moments of the eigenvalue density 
function as 



Zym = J dundun exp - l)^|Mn|^^ 



(5.11) 



where the moments, related to U by 



1 



Un ^ ^tr(C7-). 



(5.12) 
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Note that the coefficients of in ( |5.11j ) are positive for all n; it follows that Zym is 
dominated by the saddle point Un = for all n, in which the eigenvalue distribution is 
uniform. Of course, the same result is true for the eigenvalues of the holonomy matrix V. 
Thus we conclude that, independent of the torus radii, the eigenvalues of both holonomy 
matrices are unclumped for pure Yang- Mills theory on T^. 

Notice that at small A the effective mass for Un in ( |5.11J ) is very small (the mass squared 
is approximately nA^/4). As a consequence, in this regime a small attractive perturbation 
could cause the eigenvalues to clump. We will now argue that such a perturbation is 
supplied by the effective potential generated by scalar fields $/ at large but finite mass. 

5.3. Large M , noncompact limit 

We first note that when Ri is larger than any other length scale in the problem, 
the theory should behave like the noncompact thermal 1+1 dimensional theory with very 
massive scalars. This theory was analyzed in |jl^. There, it was argued that by integrating 



out all degrees of freedom except the spatially dependent holonomy matrix (in the R2 
direction) U{x), the system reduces to the model 

Using collective field theory methods, they showed that this model displays a first order 
deconfinement transition at a temperature corresponding to 

..^^Nmjfl+ofi^M)) . (5.14) 

m \ \ ln(m) / / 

Thus, our model will have a phase transition curve asymptoting to this value for large 
ri. Above the curve (corresponding to the noncompact theory at low temperatures), the 
eigenvalues of the temporal holonomy U will be unclumped, while below it, they will be 
clumped. By symmetry, there must be an additional phase boundary asymptoting to 
ri ~ ^ ln(m) at large r2, across which the eigenvalues of the spatial holonomy V clump. 

5.4. Large M , small volume 

We may now ask about the opposite limit, where Ri is so small that the theory is 
effectively one-dimensional. Assuming that the eigenvalues of V are clumped (we will see 
that this is the case for small enough Ri), the theory will behave as a one dimensional 
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theory with p massive scalars and one massless scalar (from the spatial component of the 
gauge field) as long as the Kaluza- Klein scale 1/Ri is much larger than the mass M, the 
temperature 1/ R2 and the scale -\/A of the gauge theory, yielding conditions 



ri < r2, ri ^ 1/m, ri ^ 1. (5.15) 

Here, the last relation is implied by the second for large m. 

In the language of section 4, the corresponding one-dimensional theory will have pa- 
~ i i 
rameters t = r^jr^ and m\ = mrl . From the results of section 4.4, we expect that this 

theory undergoes a phase transition at some tc{mi) for all values of mi, where tc rises 

from some value t^cF~^^^ of order one for small mi and asymptotes to 2mi/(3 In(mi)) for 

large mi. Expressing these results in terms of our two-dimensional parameters, we predict 

a phase transition at 

1 

r? 



which gives 



r2 = ^ , (5.16) 

tc{mrl ) 

. dfti^'+^^f ri < ^ 



pm 



e^"-' ^ ^ < ri < — 



where to derive the last expression we used the equation in the second line of ( [1.12|) . For 



larger values of ri, the conditions (|5.15| ) are no longer satisfied, so the theory is no longer 



well approximated by a one-dimensional theory. 

The transition line ( |5.17| ) separates a region with clumped U eigenvalues (for smaller 
r2) from a region with unclumped U eigenvalues (for larger r2). By symmetry, we will have 
an addition phase boundary separating regions of clumped and unclumped V eigenvalues, 
obtained by the replacement ri ^ r2 in ( [5.17| ). 

5.5. Large M , intermediate radius. 

At intermediate values of i?i, neither of the approximations so far apply, but the 
theory is still simple enough to treat analytically as long as the scalar masses are very 
large. In this case, we can reliably integrate out the scalars at one-loop order, leading to 
an effective action 

^ = ^ (/ rf2xtr(Ff2)) +|lndet(-D; + M2) (5.18) 
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for the gauge fields. The exphcit expression for the determinant will in general be quite 
complicated, with local terms built from the gauge field strength and its covariant deriva- 
tives, together with non-local terms involving the Wilson lines about the spatial and tem- 
poral cycles. The former will be small relative to the tree-level term as long as m ^ 1. 
The non-local terms involving the holonomy around a cycle fall off exponentially with the 
radius of the cycle in units of the inverse scalar mass, but we will see that these terms are 
important even when this exponential is very small (equal to an inverse power of m). 

For sufficiently small -Ri and R2 (we will be more explicit below), terms involving 
the spatially varying modes of the fields will be suppressed, and it is enough to consider 
the effective action assuming that the components Ai and A2 of the gauge field (and 
hence the holonomy matrices) are spatially constant. In Appendix B, we show that up to 
commutator terms, the result for the effective action in this case is 

S,ff = ^lndet{-Dl+M^) 

= -p tr(t/"y'=)tr([/-"1/-^)Mi?ii?2 

(fc,n)/(0,0) 



^ V- . V- . . ^ 2Tr^/R^k^TR^ (5.19) 



exp (-M^/R^k^l^R^ 

(fc,n)/(0,0) 



2V2^ {Rjk^ + Rln^)^ 



where in the last line we have taken the limit of large MRi and used the asymptotic form 
of the associated Bessel function Note that if MRi and MR2 are both sufficiently 

large, all terms will be small relative to the (/c, n) = (0, ±1) and (/c, n) = (±1, 0) terms. In 
addition, the commutator terms that we have not written necessarily involve at least two 
powers of U and two powers of V, so we expect them to be suppressed at least as strongly 
as the U^V'^ terms above. 

Thus, when the spatially varying modes can be ignored and both MRi and MR2 are 
large, the effective action from the scalars is well approximated by 

S,ff{U,V)=f{U)+giV), (5.20) 

where 

f{U) = -^|iv/M^e-^«Hr(t/)tr(C/t) , 

(5.21) 

g{V) = ^^^/MR'le-^^Hv{V)tv{V^) . 

V 27r Ri 



When U and V commute, this is precisely a pairwise potential between the A'' points on the 
dual torus formed by the simultaneous eigenvalues of U and V. 
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Note that / comes from the second term in ( ^.131 ) if we take the integral to run over a finite 
range Ri and ignore the non-zero modes of U{x). Both the kinetic term and the potential 
term in ( |5.13| ) will produce quadratic terms in the non-zero modes of U{x). By demanding 
that the relative coefficient for the potential term is small, we obtain an explicit condition 
for when the non-zero modes of U{x) may be ignored,^ 

rj < — ^e'"'^^ (5.22) 

Similarly, the non-zero modes of V{x) can be ignored when 

ri<— ^e"^'^!. (5.23) 



When all of our assumptions apply, the partition function becomes (using (|5.18| ) and (|5.6|) ) 

Zy^ = I DU DV J2dRe~^^'^^hR{UVU-^V-^)e-f^^h-3^^\ (5.24) 

R 

Making the change of variables U — WUW~^ for unitary W, and integrating over W 
(more details in appendix D), we find 



Zy^ = j DUDVY^e 



C2{R) ^R ^-.f(U)^-p(V) 



R R 

-1 + XR{V)XRiV^) + XRiU)XRiU^) - ^XRiU)XRiU^)XRiV)xR{V^) 

(5.25) 

For large A^, the coefficient d'j^/{d'j^ — 1) can be set to one. 

From this expression, we would like to understand the behavior of our order param- 
eters, \n{Zym)/N'^ and the expectation values for the temporal and spatial Wilson loops. 
We note first that for the purposes of computing the partition function or any observable 
depending on either U or V alone, the last term in ( ^.25|) can be ignored relative to the 
other terms as long as RiM and R2M are large enough. For, suppose we evaluate the 
partition function with an operator 0{U) inserted. Then in the last two terms of ( |5.25| ), 
the V dependent terms may be collected inside the U integral to obtain 

/ DV e-^(^) |l - ^^XRiV)xRiy^)^ = Z,\l- {^XR{V)xR{V^)),'j • (5.26) 



This condition is most appropriate in the regime ri 2> ^2 where ( 5.13| ) is vahd, but this is 
the only place we will need it in what follows. Similarly, the equation for the non-zero modes of 
V{x) is the appropriate condition for the regime r2 2> ri . 
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Here, Zg and ()g are the partition function and expectation values for the unitary matrix 
model with action giy\ But as long as the coefficient in g{y^ is small, 

r2 « -J^e^'- , (5.27) 

such a model will be dominated by the repulsive eigenvalue potential from the measure, 
and display confining behavior with a uniform eigenvalue distribution for the saddle point. 
In this case, the expectation value in ( |5.2(j| ) will vanish for large A^, since this is just the 
norm squared of the Polyakov loop in the representation K. In other words, term by term 
in the sum over K in ( p.25| ) , the contribution from the last term will be negligible relative 
to the contribution from the third term. Similarly, if we insert any operator O(l^), the 
last term will be negligible compared to the second term as long as 

ri < i./l^e'^'^^ . (5.28) 

Let us now assume that both ( |5.27D and ( p. 2 8] ) are satisfied and consider first the 
partition function (|5.25[) with no operator inserted. We then obtain 



'Z'YM - -ZM=ooZgZf + ZfZs^{V)+g + ZgZs^{U)+f (5.29) 

where Zm=oo is the partition function (|5.9[ ) for the pure Yang-Mills theory, and 5*00 is 
the action in the expression ( |5.8D for the pure Yang-Mills partition function reduced to a 
single matrix integral. Now, Zm=oo-, Zg, Zf are all finite for large (when our conditions 
are satisfied) since they correspond to confining theories, so Zym will show deconfined 
behavior if and only if either Zg^(y-^j^g or Zs^(u)+f does. 

Using the large result ( |5.10D to get an explicit expression for 5*00 and rewriting in 



terms of the modes of the eigenvalue distribution, we find 

\ ri 



^2n r2 

(5.30) 

We see that the scalar effective potential gives a negative contribution to the mass squared 
for the lowest Fourier mode of the eigenvalue distribution. Thus, ui becomes tachyonic, 
and the free energy in ( ^.30|) becomes of order A^^, when 

^ '"'v^e-'^'^^ > (l-e-^)2. (5.31) 



2tt r2 

27 



Since we assumed that the left-hand side is very small, this can only happen when rir2 ^ 1. 
Similarly, Zs^(v)+g will show deconfining behavior when 

^ ""^v^e""^^! > (l-e-^)2. (5.32) 



^27r ri 

These two curves intersect at 

4 , , w /ln(ln(m))\ , 

n = r2 = - In m l + O ) / , 5.33 

and together form part of the boundary of the region with free energy of order one and both 
U and V unclumped. Note that the curve (|5.31| ) exits the region where ( |5.22|) is satisfied 
when ri ~ m/ln(m^) and r2 has decreased to be of order 21n(m)/m. Beyond this, we 
cannot trust the functional form ( ^.31| ) of the curve, but we have already argued in section 
5.4 that there should be a phase boundary asymptoting to precisely this value, suggesting 



that our result (|5.31| ) matches on smoothly to the large Ri behavior (and similarly for 
( ^.32|) ). Noting that rir2 ^ 1 on the curve ( |5.31j ) everywhere that ( |5.22| ) is satisfied, we 
can write a more explicit expression for the phase boundary in this region. 



(5.34) 



4p i/m _ 1 
\/2iT J.2 ri 



2 



4p Vm 1 
\/27r ^2 r2 



1 

Note that the conditions (IOTD and (|]2|) are satisfied for all such values. 

To see what happens on the deconfined side of this boundary, let us now understand 
the behavior of the order parameters (|tr(t/) p/A^^) and (|tr(y) p/A^^), focusing for now on 
the first one. As we have argued above, the last term in (|5.25|) can be ignored in computing 
either of these, so by the arguments leading to ( |5.29| ), we obtain 

( J_|tr(C/)|2) = ~^M=oo^g^/(|mP)/ + ZfZs^(V)+g{\ui\'^)f + ZgZs^{U)+f{\ui\^)s^{U)+f 

N"^ -ZM=oc.ZgZf + ZfZs^{y)+g + ZgZs^[u)+f 

(5.35) 

Now, in the region of interest, at least one of Zs^(u)+f or ^Soc{v)+g behaves as exp(cA^^) 
so the first terms in both the numerator and denominator (which are at most of order one 
in A^) will be completely negligible here. In the remaining expression, 

^ ^ ^^^(^U^^2\ _ ^f^S^{V)+g{\ui\'^)f + ZgZs^^U)+f{\ui\'^)s^{U)+f 

ZfZs^(V)+g + ZgZs^{U)+f 
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even when both -^s^ ([/)+/ and Zg^(^Yj_^_g have exp(cA^^) behavior, the larger wiU com- 
pletely dominate in both numerator and denominator, so 



\ni\') 



f 







(l«ir)s^([/)+/ ^ 



^S^{V)+g > Zs^{U)+f 



(5.37) 



The model with larger magnitude for the free energy will be the one for which the \ui\ 
mode is more tachyonic, and by inspection of (|5.21|) , this will be the Soo{U) + f model for 
-Ri > i?2 and the 5*00 (^) + 9 model for R2 > Ri. 

Thus, from ( |5.37| ) and the analogous result for = (|tr("l/)p/A^^), we conclude 

that there is a phase boundary at ri = r2 starting from curve ( |5.34| ) and continuing 
towards the origin, such that the eigenvalues of U are clumped and the eigenvalues of V 
are unclumped for Ri > R2 while the opposite is true for R2 > Ri- This phase boundary 
cannot continue all the way to the origin, since we have argued in section 5.5 that both 
sets of eigenvalues are clumped for small enough ri = r2. Indeed, in ignoring the higher 
order terms in ( p.l9|) , we have assumed that RiM and R2M are large, so we can only say 
with certainty that the phase boundary exists in the region l/m^ri=r2<4 ln(m) / m. 
For smaller values of r, the phase boundary must bifurcate symmetrically, and we expect 
that the two curves thus produced connect smoothly onto the two phase boundaries which 
we argued in section 5.5 emanate from the origin. 




Figure 14: Phase diagram for Euclidean two-dimensional Yang-Mills theory 
with massive adjoint scalars on . Solid lines indicate phase boundaries for 
which we have analytic expressions. 
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5. 6. Summary of the large mass theory 

Combining all of our results, we conclude that the phase diagram for large Yang- 
Mills theory with very massive adjoint scalars compactified on a rectangular torus appears 
as in figure 14. Here, solid lines indicate regions where we have been able to determine the 
analytic form of the phase boundary. Note that this diagram displays the same qualitative 
behavior as one of the possible completions of the massless phase diagram (figure 13), so 
it seems reasonable to speculate that the phase structure is qualitatively identical for all 
values of the scalar masses. 




Figure 15: Phase diagram for 1+1 dimensional Yang-Mills theory with 
massive adjoint scalars on S^. The free energy is of order one in the shaded 
region and of order A^^ elsewhere. 



In figure 15, we present an alternate version of the diagram, with axes labeled by the 
size of the spatial circle in units of the coupling (ri = R\/\) and the temperature in units 
of the spatial radius {TR = ri/r2). Analogous units were used in our analysis [p|,p!5[] of 
3+1 dimensional Yang-Mills theory on a sphere, where our analysis suggested a single first 
order transition extending from large volume to a finite value of TR in the zero volume 
limit. In the present case, it is interesting to note that: 

1) The solid line corresponding to deconfinement as measured by the temporal Wilson 
loop expectation value is not smooth. 

2) While the deconfinement transition does extend all the way to zero volume as for 
the topologically trivial case, the transition temperature does not approach a finite value 
in units of the inverse spatial radius. 
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3) The large volume regions for both the confined and deconfined regions are separated 
by phase transitions from the small volume regions. 

4) The two order parameters for confinement, l/N"^ times the free energy and the 
expectation value of the Polyakov loop, do not have the same behavior at small and large 
volumes in this example. More specifically, in the "confined" region below the solid line, 
the free energy is of order one at large volume, but it becomes of order A^'^ to the left of 
the dotted line. 



6. Maximally Supersymmetric Quantum Mechanics 

In the rest of this paper we will study one and two dimensional maximally supersym- 
metric gauge theories on circles and tori. As is well known, Yang-Mills theories with 16 
supercharges display rather different dynamics from their bosonic counterparts; in partic- 
ular the scalar effective potential in these theories has fiat directions, and the spectrum of 
these theories is not gapped. As a consequence we expect (see section 2) that the phase 
diagrams of these theories (which we will obtain) will differ in significant aspects from the 
equivalent diagrams for the purely bosonic theories of the previous sections. 

Partition functions of supersymmetric theories on tori are strongly dependent on the 
boundary conditions of the fermions around cycles of the tori. When all fermions are 
periodic, the path integral is localized on supersymmetric configurations (and computes 
a Witten Index [|16[)- These path integrals depend weakly (if at all) on parameters, and 
do not undergo phase transitions. In this paper we will focus on the dynamically more 
interesting path integrals with anti-periodic boundary conditions for fermions; these are 
relevant in particular for studying finite temperature. 

In this section we consider the maximally supersymmetric SU{N) Yang-Mills theory 
in one dimension, in the 't Hooft large limitEl. The Euclidean Lagrangian is given by 

S=y^jdt tr{^Do<5^Do<5^ + '-^^ Dq^ - ^ + ^) , (6.1) 

where z, j = 1, ■ ■ ■ , 9 and \1/ is a real 16 component spinor. The time direction is identified 
with a period /3 = 1/T; the bosonic fields are periodic in this direction and the fermions 



It is important to distinguish this from the M(atrix) theory large A'' Umit |17| , which is 
conjectured to describe the hght-cone quantization of a flat-space gravitational theory (M-theory) 
and has a different thermodynamic behavior. 
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are anti-periodic. The theory may be characterized by the single dimensionless parameter 
t = T/A^^, in terms of which the effective dimensionless 't Hooft coupling is l/P. 

As the fermions are anti-periodic about the circle, all fermion modes acquire large 
masses as the circle shrinks to zero size, so in the small radius or large t limit, our the- 
ory behaves identically to the bosonic theory considered in section 3; in particular, the 
eigenvalues of the holonomy matrix are clumped on the circle. 

In the opposite, small t, limit the quantum mechanics is strongly coupled; but in this 
case, we can use a dual description to understand the dynamics. 

We recall that the 't Hooft limit of massless maximally supersymmetric 0-1-1 dimen- 



sional Yang-Mills theory is believed [lT8| , p!9|] to be dual to type IIA string theory in the 
near horizon limit of the extremal black DO-brane solution of type IIA string theory. This 
background has a reliable description in terms of type IIA supergravity (small curvatures 
in string units and weak string coupling) for a range of the radial coordinate corresponding 
to energy scales E <^ X\^^ = (gfN)^ in the quantum mechanics.0 

Similarly, the finite temperature quantum mechanics in the 't Hooft limit should be 
dual to a near-extremal version of this geometry. This dual geometry contains a horizon, 
and type IIA supergravity is a valid description near the horizon as long as t ^ 1. The 
presence of a horizon in the Lorentzian geometry implies a contractible thermal circle in 
the associated Euclidean geometry. Thus, the expectation value of the (traced) Wilson 
line around the circle in the gauge theory, which is mapped in the string theory (after 
adding to it some dependence on scalar fields) to the action of a string worldsheet ending 
on the thermal circle ||2^,^,22] is non-zero. This implies that the eigenvalue distribution 
is non-uniform. 

Therefore, for the massless supersymmetric theory, the eigenvalues of the Wilson line 
operator seem to be clumped (or at least non-uniformly distributed) at both large and 
small values of t, so the phase transition that we found in the bosonic theory appears to 
be absent here. 

This disparity of behavior between the bosonic and supersymmetric theories at low 
temperatures should not come as a surprise; indeed, their zero temperature dynamics 
is completely different. In the purely bosonic theory, the classical flat directions in the 
potential are completely removed by quantum fluctuations which result in a linear potential 



For large but finite N, the type IIA supergravity description is also no longer valid at low 
energies, when E <C Qi N ' , but this scale goes to zero in the 't Hooft large A'' limit. 
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between the spatial eigenvalues. The resulting spectrum should be discrete. On the other 
hand, for the supersymmetric theory, the quantum fluctuations generate only the well- 
known x'^/lxl^ potential between eigenvalues, and the spectrum is continuous. 

As for the bosonic theory, we can consider more generally a massive deformation. 
The simplest such deformation has equal masses for all scalars and fermions (of course, 
this breaks supersymmetry). In the limit of large mass, all fermion and scalar modes are 
weakly coupled, so as for the bosonic theory, we may integrate out the matter at one 
loop to reduce the partition function to a unitary matrix model. The result is similar to 
( [4.6|) , but with p = 17 (assuming that all fields have the same mass). Thus, for large 
enough masses, we again have a phase transition from clumped eigenvalues at large t to a 
uniform distribution for small t, with tc = m/ln(17) + (9(l/m^). Since we have a phase 
transition for t ^ m/ln(17) at large m and no phase transition at m = 0, it must be 
that the phase transition line intersects the t = axis either at m = or at some small 
non-zero value of m. There is certainly a qualitative change in the theory as soon as a 
mass is turned on, since the zero temperature potential goes from being flat asymptotically 
(with a corresponding continuous spectrum) to asymptotically harmonic (with a discrete 
spectrum). While this motivates a possible phase transition, it is not clear how to directly 
relate this information to the behavior of our order parameter. 

7. Maximally Supersymmetric Yang-Mills Theory on 

In this section, we would like to study the thermodynamics of maximally supersym- 
metric l-fl dimensional Yang-Mills theory on a circle, corresponding to the Euclidean 
theory on T^. In this case, the fermions are anti-periodic about the thermal circle, but we 
have a choice between anti-periodic and periodic boundary conditions for the spatial circle. 
In fact, we will see that both possibilities are included if we make the natural generalization 
to non-rectangular tori. 

7.1. General tori: classification and fundamental regions 

In this subsection we will classify and describe the inequivalent tori on which max- 
imally supersymmetric Yang-Mills theory can be compactified. We specify our torus by 
two identifications: 
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(where r = ti + zt2 is the modular parameter), and further specify that all fermions are 
anti-periodic around both of these cycles. Below we will denote the holonomy along the 
cycle z ^ z + L as V and the holonomy along z z + Lt as U. The phase diagram we 
wish to determine is a function of three real dimensionless variables: the complex variable 
T and the dimensionless coupling constant A = gyM^^-i where A = L^T2 is the area of 
the torus. 




Figure 16: 

fermions. 



Fundamental region of the modular group for a theory with 



Bosonic Yang-Mills theories on tori whose r parameters are related by an SL{2, %) 
transformation (r (ar + b)/{cT + d) with integers a,h,c,d satisfying ad — he = 1) are 
identical at equal values of A. It is thus sufficient to study such theories in the familiar 
SL{2, fundamental domain, denoted by the region I in figure 16. 

The supersymmetric Yang- Mills theory we will study in this section includes fermions; 
the partition function for this theory will be identical (at equal values of A) only on those 
tori that are related by the subgroup of SL{2,%} that preserves our fermion boundary 
conditions. This group is generated by r — > — l/r and r — > r -|- 2. The fundamental 
domain for this subgroup is the union of the regions I, II and IV in figure 16. In this 
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section we will study the Yang- Mills partition function on tori with modular parameters 
lying in this fundamental domain. 

To end this subsection we comment on the physical interpretation of tori in regions II 
and IV. Tori on the two vertical lines that border the diagram in figure 16 have modular 
parameters of the form r = ±1 + iT2 (with real T2). These tori are best thought of as 
rectangular tori with r' = iT2, with periodic boundary conditions on fermions along the 
r' cycle (the fermion boundary conditions remain anti-periodic along the '1' cycle). In 
general, a torus in the region I lb/ a with modular parameter r may be thought of as a 
torus in the region I, with r' = r ± 1 and periodic boundary conditions along the r' cycle. 
Similarly, a torus in the region IVb/a may be reinterpreted as a torus in the region I with 
t' = — 1/(t±1) and periodic boundary conditions along the '1' cycle. 

7.2. Analysis at small A 

As discussed in previous sections, at small A the Kaluza Klein modes of the field 
theory on may reliably be integrated out at one loop, provided that A/M^ is small, 
where M is the mass of the lightest KK mode in units in which the area of the torus is 
one. When this integrating out is legitimate, the resulting determinant (see Appendix A) 
ensures that the eigenvalues for both holonomies U and V clump. 

Expanding all fields (take a scalar field (/){z) as an example) on the torus in Fourier 
modes subject to the periodicity (j){z) = (j){z + L) — (j){z + Lt) gives 



9m,n exp 
m,n=— 00 



27ri 

- — (T2mRe(2;) + (n — Tim)Im(^)) 

Lt2 



(7.2) 



leading to = ^ (m^|Tp + m? — 2mnTi) (in units of l/L). For bosonic fields m and n 
are integers, while for fermions m+1/2 and n+1/2 are integers. 

In regions I and II the lightest non-zero-mode for scalars is the (m = 0, n = ±1) mode 
with = 47r^/T2. So, we can integrate out all the KK modes at one loop reliably in 
regions I and II if 

At| < 1. (7.3) 

For larger values of T2 all modes with m 7^ (the boundary conditions ensure that this 
includes all the fermions) are still very heavy and decouple, so we can study the theory 
by including only the m = modes. The resulting quantum mechanical theory is simply 
the bosonic quantum mechanics of section 4, with R = Lt2 and coupling Ai = A/L, giving 
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an effective dimensionless coupling Ar|. It follows from the analysis of section 4 that the 
system undergoes a phase transition at T2 = a/^/x, where a is a number of order unity, 
determined in the analysis of section 4. At larger values of T2 one of the ^at symmetries 
of the system is restored, and (tr([/)) = 0. 

In region IVa the lightest mode is (m = n = ±1), with mass 



r2 



(7.4) 



(the IVb region is related to this by (ri — ri, m 
coupled and can be integrated out provided that 



-m)). The KK modes are weakly 



((i-n)2 + r|)2 



< 1. 



(7.5) 



When (|7.5| ) is not obeyed (but still A ^ 1), all modes with n ^ m become very heavy and 
decouple. The resulting effective quantum mechanical theory of the light modes includes 
fermions, and is, in fact, precisely the one dimensional supersymmetric quantum mechanics 
of section 6. We have argued that it is unlikely that this system undergoes a phase 
transition. It seems that our system is in the (tr(t/)) 7^ 0, {tr{V)) 7^ phase throughout 
region IV. 




Figure 17: The phase transition line at T2 ~ 1/V A and its images, for a 
specific value of A. 
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Mapping regions I and II by the (allowed) modular transformations we find an infinite 
number of images of the phase transition line T2 = a/\/A, which are drawn in figure 17. 
For example, in region III this line is mapped by the r ^ — 1/r transformation to the line 
— \/A/2a)^ + t'i = X/Aa^. Above this line both holonomies are clumped, {ix{U)) 7^ 
and (tr(y)) 7^ 0, while below the line (in region III) (tr(y)) = 0. 

To summarize, a torus is characterized by its modular parameter r. If we scale 
the 't Hooft coupling constant appropriately then the partition function of the theory 
is a modular invariant function on the r plane. Thus, we only need to study it in the 
fundamental region (which is different from the standard one because of the fermions) . In 
this region we find a single phase transition line where the eigenvalues of the Wilson loop 
along the longer fundamental cycle clump. This line has an infinite number of images on 
the r plane through modular transformations. Above we gave the results for the case of 
two anti-periodic cycles for the fermions, but the case of one periodic cycle (corresponding, 
in particular, to the finite temperature theory) is simply related to this by r — > r + 1. 

7. 3. Strong coupling from AdS/ CFT 

When the effective coupling A is large, the theory cannot be analyzed using perturba- 
tive techniques, but we can use the AdS/CFT correspondence to study the phase structure 
of the system. 

The dual of the 1 + 1 dimensional maximally supersymmetric Yang-Mills theory in 
the 't Hooft large limit is the near- horizon geometry of the Dl-brane metric of type 
IIB string theory |jl^. If we put the gauge theory on a spatial circle of circumference 
L then the dual theory is obtained by a periodic identification of the spatial coordinate. 
If we want to study the theory at temperature T we have to consider the near extremal 
Dl-brane solution. The relevant Euclidean string frame metric and dilaton are: 



ds^ = a' 



^3 ul\ df" d6'^ 



1-^ -^ + ^==-—^ + 



VdJ7\ u^J yrf^(27r)2 ^3(i_< 



^,. du' + u-'^/d^dnj 



(7.6) 

where A' = AL^, di = 2^n^, = i^^TL^V, we have to identify t' ~ t' + 1/T to avoid 
the conical singularity at u = uq, and the periodicity of 9' is 2n. There is also a 3- form 
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field which we do not exphcitly write down. Defining p = TL and changing coordinates to 
t = t'T and 6 = we can write the metric in the form: 



a 



Vd^P^ 



Ur 



dr + d9' 



+ 



Vd^ 



du^ + u-^ ^d^ dO^j 



(7. 



with Mq 



1671- 



5/2 



3 -pyX'. Now, defining a complex coordinate ^ = t + z^, we have the 
identifications z ^ z -\- 1 ^ z -\- ip, and the conformal boundary of this geometry is a 
rectangular . The parameter p appears both in the identifications and explicitly in 
( [7.8|) , but in fact we can get rid of it in the metric by noting that the form of (|7.8| ) is 
invariant under a rescaling u au, uq au, X' ^ a^X' , p ^ ap. If we use this rescaling 
with a = we can rewrite ( [7 .81 ) as 



ds' 



a 



diX L 



1 - 



dr + dO' 



+ 



diX 



1 _ i!o 



du^ + u VdiAdO^ 



(7.9) 



with A = A'/p^ and Uq 



167r 



5/2 



g V Since now p appears only in the identifications, it 
is clear that if we change the identifications toz~z + l~z + r with r = ri + zr2 an 
arbitrary complex number, we still have a solution to IIB supergravity (with the same uq) 
without any singularities, whose asymptotic boundary is a torus of modular parameter r. 

If we want to compute the thermal partition function of the field theory living on the 
boundary we have to consider all supergravity (or more generally string theory) solutions 
Xi with the appropriate behavior at infinity. The partition function of the dual gauge 
theory in the large limit is: 



Zgauge = ^ exp( -/(X^ ) ) , 



(7.10) 



where / is the Euclidean action of the supergravity (string theory) solution. 

Let's say that we want to compute the partition function of supersymmetric Yang- 
Mills (SYM) theory on a torus of modular parameter r/. It is clear that solutions of the 
form ( |7.9| ) contribute to the partition function only if the asymptotic of the supergravity 
solution can be conformally mapped to the on which the gauge theory lives. This is 
possible if and only if the modular parameter of the torus of the gravity Tg is related to Tf 
by some modular transformation. So we have to start from the supergravity solution with 
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Tg = Tf, consider all modular transformations of Tg, and then sum the exponentials of the 
Euclidean action of the corresponding supergravity solutions P1,B3 . 



In general all solutions will contribute to the partition function. But we notice that the 
Euclidean action scales as A^^, so in the large limit we get the dominant contribution from 
the solution with the smallest Euclidean action. This allows for sharp phase transitions in 
the N ^ oo limit, if the solution with the minimum action changes discontinuously as we 
vary the parameters of the theory. 

If we try to compute the Euclidean action of any of the solutions we find that it 
diverges as we integrate over all space. However we are not interested in the value of 
the action but in comparing the action of different solutions. One way to compare the 
actions is to match the solutions (metric and other fields) at some large cutoff value of 
the radius, calculate the action difference and then send the value of the cutoff to infinity. 
The divergent terms cancel once the geometries are properly matched, and one obtains a 
finite result. Alternatively one can add counter-terms to cancel the divergences, leading 
to equivalent results. 

Before we proceed we must clarify two subtle points: 

(a) Not all solutions with different Tg are different. For example, it is easy to see that 
Tg and Tg ± 1 describe the same supergravity solution. The asymptotic torus of the 
supergravity solutions has a cycle that is special, namely the (1,0) cycle (labeled by 
t) which is contractible in the interior. It is fairly easy to see that all inequivalent 
geometries can be characterized by specifying the cycle (p, q) of the gauge theory torus 
onto which the (1,0) cycle of ( [7.9|) is mapped when we conformally match the two 
tori. The two integers p and q have to be relatively prime. 

(b) We have to be careful about the boundary conditions of the fermions and the spin 
structures of the supergravity solutions. If there is a circle factor in the asymptotic 
geometry then we can choose either periodic or anti-periodic boundary conditions for 
the fermions along this circle if it is not contractible in the interior. However, if it 
is contractible then only anti-periodic boundary conditions are allowed. In our case 
the (1,0) cycle of the gravity torus is contractible, so it has to be mapped to an anti- 
periodic cycle of the gauge torus. Since in our gauge theory analysis above we chose 
anti-periodic boundary conditions for the fermions on the (1,0) and (0,1) cycles of 
the gauge theory, we conclude that the acceptable geometries of (a) are those with 
p + q odd. 
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Let us now analyze which solutions we have to include for a specific gauge theory on 
a torus with parameter r whose partition function includes a contribution from ( [7.9| ). Up 
to a change of coordinates all these solutions should have the same asymptotics as ([7.9|), 
and they should involve identifications by some modular parameter r' which is related 
by an SL{2, %) transformation to r. It is easy to verify that for this to be the case the 
parameters of the two metrics have to be related by AV2 = Ar2, and the u coordinates 
of the two metrics are related by u' ^fi^ = Uy/r^. Note that the relation between the A's 
is precisely the one we expect from the field theory point of view, for the two theories to 
have the coupling constant times area. 

A straightforward computation of the (regularized) Euclidean action of the IIB solu- 
tion ( |7.9| ) gives: 

I = ^-^^T2 / duu^, (7.11) 
64A27r9 Juo 

where R is some very large radial position and is the volume of the unit 7-sphere. Using 
the results of the previous paragraph we see that if we compute I{Xi) — /(X2) for two 
solutions with r and r' related by a modular transformation, the divergent parts cancel 
and we find: 

I(X,) - I(X,) = + ^) , (7.12) 



9\/AT2?r3/2 

where we used Uq = — V A. So, using the invariance of Ar2, we conclude that the 
solution with the minimum action is the one which has the maximum value of ti (and is 
consistent with p + q = odd) . 

Let us now see how this works in different regions of the r/ plane. In regions I and 
II the geometry that dominates is Tg = Tf. It maps the (1,0) cycle of the gravity torus 
to the (1,0) cycle of the gauge torus. Any other solution with Tg related by a modular 
transformation has a smaller T2 so it has a bigger action. Since the (1, 0) cycle is the only 
contractible cycle in the gravity solution, the Wilson loop around it will generically be 
nonzero (tr(t/)) 7^ 0, while all Wilson loops around any other cycle will be zero. Though 
it takes a little more work to see it, region IV is also dominated by the same geometry^. 

So, in summary, a single saddle point ( |7.8D dominates the thermodynamics of strongly 
coupled Yang-Mills theory on a torus in the fundamental region (the union of /, //, IV). 



In region IVa there is the geometry with the largest T2 is the Tg = — .However this would 
map the contractible (1, 0) cycle to (1, ±1) which is a periodic cycle, so this is not an acceptable 
solution. 



40 



The (1,0) cycle is contractible on this solution, so (tr(y)) = and (tr(t/)) in this 
phase. For other regions we can find the dominant solution by mapping them into the 
fundamental region. 

l.Jt.. Putting it together 

Let us summarize our understanding of the phase structure of this system. We begin 
by discussing tori with r = it where t is real. At weak coupling the system is in the 
(tr(t/)) ^ 0, (tr(F)) = phase at large t. At t = a/^/x the system undergoes a phase 
transition; at smaller values of t the system is in the (tr(t/)) ^ 0, (tr(F)) ^ phase. At 
t = 's/x/a the system undergoes another phase transition, to the (tr(?7)) = 0, (tr(F)) ^ 
phase. On the other hand, at strong coupling the analysis of the previous subsection shows 
that the system undergoes exactly one phase transition at t = 1. When t > 1, (tr(^7)) ^ 
and (tr(F)) = 0. On the other hand, when t < 1 (tr(y)) ^ and (tr(t/)) = 0. 
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Figure 18: Conjectured phase diagram for SYM on a rectangular torus 
with anti-periodic boundary conditions on both cycles. 

The phase diagram in figure 18 summarizes this behavior, and provides the simplest 
possible interpolation between these two limits. 
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(antiperiodic) 



^^^.^ large X 

ti-(U)*0 
tr(V)=0 



tr(U)*0 /- ^ 

tr(V)*0 / 

/ ^ ^ ^ - ~ ^ ^ small X t 



r I (periodic) 

Figure 19: Phase diagram for SYM on a torus witli one periodic and one 
anti-periodic boundary condition for the fermions. The large A transition 
here is seen in the dual IIA description as a Gregory-Laflamme transition 0] . 

Next, recall that a torus with Re(r) = 1 may equally well be regarded as a rectangular 
torus with r' = r — 1 but with periodic boundary conditions along the r axis. Following 
previous discussion in we expect a phase transition associated only with the V holonomy, 
seen as a Gregory-Laflamme transition in the type IIA string description which is obtained 
by performing a T-duality on the periodic torus direction. Hence, we expect the phase 
diagram of figure 19 (or equivalently figure 1 in From the analysis in |^] we expect the 
phase transition to occur at t ~ 1/A for large A. Note that this phase transition involves 
a new solution (localized in the T-dual direction) which is not included in the discussion 
of the previous subsection. 

Using this information we conjecture the phase structure for a general value of r. In 
figure 20 we show the two phase boundaries we expect for clumping of the V holonomies 
in the fundamental domain as a function of r and A. For large A, the relevant boundary 
(shown in red) coincides with the boundary of the fundamental domain and corresponds 
also to anti-clumping of the U holonomy. For small A, the relevant boundary (shown 
in blue) involves only the transition of the V holonomy. As discussed earlier, modular 
transformations of this latter surface give other phase surfaces associated with transitions 
of different holonomies; for example, the clumping of the U holonomy is given by r — I/t. 
The Re(r) = slice reproduces the previous figure 18 when converting to ri, r2 coordinates 
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Figure 20: Conjectured phase boundaries for tlie V liolonomy siiown in the 
fundamental domain for general values of A. The phase transition associated 
with the blue surface involves only the V holonomy, whilst the red surface 
corresponds also to a transition of the U holonomy. 

and including this U transition. Likewise, the Re(r) = 1 slice reproduces figure 19. Note 
that we expect that the two phase transition surfaces (red and blue) do not join smoothly 
for any Re(r), and their intersection moves to larger A as we move nearer Re(r) = 1, 
eventually leaving only the two phases with clumped U holonomy that are seen in [^. 

By comparing figures 18 and 19 with the phase diagrams 13 and 14 of section 5, it 
is clear that for either choice of spatial boundary conditions, the supersymmetric theory 
has qualitatively different behavior from the bosonic theory. Curiously, the bosonic the- 
ory and the supersymmetric theory with anti-periodic-anti-periodic, anti-periodic-periodic, 
and periodic-periodic boundary conditions, have four, three, two, and one distinct phases, 
respectively. 
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Appendix A. One-Loop Effective Potentials for Zero Modes 

In this appendix, we shall derive the one-loop effective action for the zero modes 
of SU{N) gauge theories with adjoint scalar fields in various dimensions, by integrating 
out all massive modes. We begin in §A.l with an analysis of the effective action for the 
diagonal modes in the 0-dimensional matrix integral and its implications for the matrix 
integral discussed in §3. The same analysis applies also to the KK zero modes of higher 
dimensional theories; in §A.2 we add the contribution of the higher KK modes as well, and 
in §A.3 we consider the one-loop contribution for supersymmetric Yang-Mills theories. 

A.l. The zero mode integral 

In this subsection we investigate the effective action for the matrix integral with p 
massless scalar fields discussed in section 3, 



where a,/? = 1, • • • ,p. The same integral (except that some of the '^"'s become periodic) 
arises for the zero modes of a d-dimensional gauge theory on with pa scalar fields and 
coupling constant A^; in this context p = pd + d since additional scalar fields arise from 
the gauge field, and Aq is A^ divided by the volume of the T^. 




(A.l) 
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A naive extrapolation of the formula (|3.4|) , giving the width of the eigenvalue distri- 
bution for large masses, to small masses would suggest that the eigenvalue distribution is 
infinitely spread out at m = 0. This estimate is too crude since it ignores the effect of 
the quartic restoring term in (|A.1|) . Taking this term into account, a simple dimensional 

analysis suggests that the characteristic width scale for the eigenvalue distribution at the 

1 

saddle point for m = is given by a = KXq, where K is independent of and Xq. The 
N dependence of this estimate follows from 't Hooft scaling, while the Ao-dependence is 
deduced from a simple change of variables in (|A.1|) . However, this scaling argument can- 



not rule out the possibility that K = oo, a danger that seems real given the fact that the 
quartic restoring term in ( |A.lj ) vanishes along a noncompact 'moduli space' along which 
all the matrices commute (and so may simultaneously be diagonalized). Consequently, 
the question of whether or not the eigenvalues in ( |A.1|) clump at m = is a dynamical 
issue 0. One way to analyze this is by an analysis of the quantum effective action on 



the classical moduli space [^,|2^,^ that is obtained by integrating out the off diagonal 
modes. The generic effective mass of an off-diagonal mode is of order a, the width of the 
eigenvalue distribution. Consequently, as in the analysis of the massive theory in section 3, 
the effective coupling for the off diagonal modes is ~ A/a^; these modes may accurately be 
integrated out at one-loop provided that this coupling is small. This one-loop computation 
was performed in and is reproduced, for the convenience of the reader, in the following 
paragraphs. 

Let us expand around a configuration in which all matrices i/'" are diagonal with eigen- 
values 'i/'f (z = 1, ■ ■ ■ , A), and integrate out the off-diagonal components of the matrices in 



( |A.1| ) at 1-loop. Expanding the action ( |A.1| ) to quadratic order in the off-diagonal modes. 



we find 

^rot:3: = ^ E E ^S^f^ - (^^").. (^^')., ^sv'fj , (a.2) 

where we defined (A'i/'")ij = — '0". Naively, integrating out the (ij)'th off-diagonal 
components of all the matrices simply yields det(M)~^, where M is the matrix 

M^P = S'^^J2(Ar)l - {A,l;^hj{A,l;%. (A.3) 

7 



16 



See [25 1 for a rigorous proof of the convergence of this integral, and [26| and references therein 



for a nice review and more details. 



45 



However, the eigenvalues of the matrix M are easily seen to be with degen- 

eracy p — 1 and zero with degeneracy one. The appearance of the zero eigenvalue is a 
result of the fact that we have neglected the existence of flat directions in the space of 
off-diagonal components, due to the remaining SU{N) gauge symmetry of ( |A.1[ ), which 
can be used to rotate the diagonal components into off-diagonal components. 

To correct for this, we eliminate these flat directions by using the SU{N) symmetry 
to diagonalize one of the say i/^^, exactly. Performing this diagonalization by inserting 
a gauge-fixing 5-f unction introduces the Fadeev-Popov determinant factor 



det ~ n (^^')^, ~ n(^^') 

i^j i<j 



(A.4) 



where we have evaluated the determinant at 1-loop. 

Note that this factor is nothing more than the Vandermonde determinant appearing 
in the change of variables from the matrix i/j^ to its eigenvalues. 

Now there are no longer any off-diagonal components for in (|A.2D . This has the 
effect of removing the first row and column from the matrix M in (|A.3| ) . The eigenvalues of 
the resulting matrix are [^ji'^''P'^)ij] with degeneracy p — 2 and [Atp^)'^^ with degeneracy 
one. Having computed the relevant determinant, we may now write the contribution to the 
path integral from the gauge-fixing and from integrating out the off-diagonal components 
of the zero modes as 



p-2 



(A.5) 



The Fadeev-Popov factor in the numerator cancels the contribution from the eigenvalues 
depending only on '0^, leaving us with an expression which (as expected) is symmetric 
under permutations of the "(/^"'s, 

2-p 



n 

i<j 



(A.6) 



Thus, we find that the one-loop contribution to the matrix integral may be expressed 
as an integral over the eigenvalues of the matrices, 

1 



n.<JE,(V'7-^j)' 



(p-2) 



valid when a ^ Aq. 



(A.7) 
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In the large limit it is natural to rewrite this in terms of an eigenvalue distribution func- 
tion p{i{^°'), and to evaluate the integral by saddle points (since the action is proportional 
to N^). 

As all pairwise eigenvalue forces are attractive, the minimum of the action is attained 
when all eigenvalues sit at a point. Unfortunately this saddle point lies outside the domain 
of validity of ( |A.7|) , since it has a = 0. Thus, the one-loop analysis does not suffice to 
analyze the matrix integral, though it does reliably show that the eigenvalues attract each 
other at large distances. A more precise analysis shows that on the saddle point of ([A.l| ) 
the eigenvalues have a sharply localized distribution with the scale a ~ Aq , as suggested 
by the scaling arguments presented earlier in this section. This is verified by the Monte 
Carlo simulations shown in section 3, and was originally shown by similar methods in ||29|| . 

In summary, the behavior of ( |A.1D is qualitatively unaffected as the effective coupling 
constant Xo/m^ is varied from zero to infinity. The eigenvalue distribution of the matrix 
jj-{J2ai'^°')'^) (|A.1| ) is governed by a saddle point that is always strictly localized; the 
localization length is given by = ^^/(^): where /(O) = 1 and f{x) oc 1/y/x at large 

X. 



A. 2. Integrating out KK modes at 1-loop 

We now consider a (i-dimensional gauge theory compactified on a torus T'^ with circles 
of circumference {fJ' = 1, ■ ■ ■ , c?), and compute the effect of integrating out the KK modes 
at 1-loop on the effective action for the diagonal zero modes which we computed in the 
previous subsection. Consider the d-dimensional bosonic action0 



S 



d X -— tr 
4A 



I I, J 



(A.8) 



with p — d scalar fields, such that the zero mode action is given by ( A.l ), where t/j'^ 
(n = 1, ■ ■ ■ , d) are the zero modes of A^, and tp'^^^ (I = - - iP — d) are the zero modes 
of (j)^ . We expand the fields A^, (f)^ in KK modes A^^^^^y ^j, 't'lmx} ij^ with mode numbers 



If there are also fermionic fields, they decouple from the bosons at 1-loop so we can consider 
them separately as we will do below. 
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mx G ^ around the A'th circle. Expanding the action to quadratic order, we obtain 



N 



{mt^} i<j 



+E 



I A |2 

\^u,{mx},ij I 



\<Pi 



- E (i^r), - ^) 



A A* 

^M,{"^A},«J^!^,{mA},ji 



fi,i 



(A.9) 

Let us now perform a naive first attempt at integrating out the KK modes at one-loop. 



Defining a vector D 



ik,{mx} 



as 



27rm 



R, 



(A.IO) 



we see that performing the quadratic integral over KK modes in ( |A.9| ) yields a factor of 
det(M')-\ with 

M' = (D^D) I - DD'^, (A.ll) 

similar to what we found in the zero dimensional case. This determinant is easily evaluated, 
as the eigenvalues of M' are again D^D with degeneracy p — 1 and zero with degeneracy 
one. The appearance of a zero eigenvalue is not surprising. Indeed, it is expected as we 
have yet to fix the gauge. A convenient set of gauge-fixing constraints to adopt is the 
following : 

diAi =0, 



82 / dxi A2 =0, 



^3 / dxi dx2 As =0, 



(A.12) 



dd J dxi. . .dxd-iAd =' 
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0. 



For each constraint, we must also insert an appropriate Fadeev-Popov determinant factor. 
For the generic constraint 

dn+i j dxidx2...dxnAn+i = Q, (A.13) 
the corresponding determinant takes the form 

(n) (n)' 

det {dn) det {d^+i - i [A^+i, *]) , (A. 14) 

where the superscript (n) is used to indicate that the determinant is taken over ah modes of 
the gauge field that are constant in xi, . . . , x^, while the prime in the second determinant 
indicates that it also includes only non-zero modes in Xn+i- We neglect the first, constant, 
determinant and evaluate the second one at one-loop, finding 

m{m, - ^) ^ nn ((A*--). - (a-) 

where the product over is suitably constrained. Summarizing, we find that our gauge- 
fixing procedure introduces a factor of (^{Ai/j'^) -j — 2Tvm'^ / R^^ for each mode of the gauge 
field which has = ■ ■ ■ = mJ^~^ = and rn^ 7^ 0. 

Returning now to the task of integrating out the KK modes in ( A. 9 ), we write the 
gauge fixing constraints (|A.12|) as 



^r,m2,m3,...,md ^ 
^0,r,m3,...,md 



(A.16) 



A'^ — n 

^0,0,...,0,r — 



where r 7^ and the other are arbitrary. It is now easy to see that, for every choice 
of mode numbers {m\}, exactly one component of the gauge field, say Aj^ s^j^^^, is elim- 
inated by these constraints. When the corresponding row and column are removed from 
M', the eigenvalues become D^D with degeneracy p — 2 and ^(A'0^)j^ — 2TTm'^ / with 

degeneracy one. Each factor of ^(A'i/''')^^ — 2Txmy / R^ that arises for a gauge field com- 
ponent eliminated by the gauge-fixing conditions, though, serves to cancel a corresponding 
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factor from the Fadeev-Popov determinants (|A.15| ), so that the final result obtained from 
gauge-fixing and integrating out KK modes at 1-loop is given by exp{—Seff,bos), where 



Seff,bos = ip-2) J2 



(A.17) 



We have included also the terms (|A.6|) obtained from integrating out the off diagonal zero- 
modes. Like ( [A.Tp , the result corresponds to a pairwise logarithmic effective potential 
between the eigenvalues, but now these live on the dual torus in the gauge field directions, 
so we have additional interactions between the eigenvalues and the infinite set of image 
eigenvalues. These image interactions (corresponding to the sum over m) ensure that the 
full result is periodic. This sum may be evaluated explicitly using the results of appendix 
B, and we find 



Seff,bos = -2{l[R,)ip - 2) J]</^(27r)-^/^ 



E 



d/A 



(A,18) 



where 9j- is defined as 



(A. 19) 



A. 3. Generalization to supersymmetric field theories 

It is a simple matter to include also fermionic fields, if they exist. Assuming that the 
Yukawa couplings to the scalars are the same as the gauge couplings, as in supersymmetric 
gauge theories, and that the fermions are periodic around {d — 1) of the circles but anti- 
periodic around the d'th circle, the appropriate determinant is easily computed and found 
to yield 



nn 



d-1 



p-2 



■(A.20) 

As with (|A.17|) , this corresponds to a pairwise potential between each eigenvalue and an 
infinite series of "images" of the remaining eigenvalues, except that the potential has the 



50 



opposite sign and the coordinates of the images are shifted by a half period in any direction 
for which the fermions are anti-periodic. 

Combining ( |A.(3| ), ( |A.17D , and ( |A.2UD , as appropriate for a supersymmetric gauge 



theory on a torus, we arrive at the fuU result : 



m'^ i<,j 



2 



ij 



E';; ((AV'^ - ^) +((A^%-^) + (AV'^+^)?, 

(A.21) 

Using the results of Appendix B we can again evaluate the product over mode numbers 
and obtain an effective action for the eigenvalue separations. The contribution from the 
fermionic sector is identical to the bosonic result (|A.18|) , with the exception of an additional 



factor of (— in the sum. As a result, the fermions serve to eliminate half of the terms 
in the kd sum. Our final result for Seff thus becomes: 



^e// = -4(n^M)b-2)E< (2 



i<j 



fc^=-oo fcd odd klRl 



(A.22) 



A. 4- Analysis of the effective potentials for d = 2 

We will now focus on the case d = 2, and argue that for both the bosonic and the 
supersymmetric cases, the eigenvalues of the gauge field will be clumped in all directions 
whenever the one-loop effective potentials ( [A.17P or ( |A.21| ) are reliable. 



The contribution of the zero modes (corresponding to the m = terms in the bosonic 
contribution to the potential) was discussed in §A.l, and we concluded that this leads the 
eigenvalues to clump at a scale 

a ~ (A.23) 

where Aq = A/(i?ii?2)- The additional terms in the effective potential correspond to the 
interactions with image eigenvalues. These interactions will be negligible for the clumped 
configuration as long as the scale (|A.23| ) is much smaller than either 1/Ri or I/R2 (i.e. 
when the distance to the nearest image charge is large compared with the distances between 
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eigenvalues in the clump). Il^ Since the KK-mode terms are important for configurations 
that are not highly clumped, we may worry that the full effective potentials ( |A.l^ ) or 
( |A.22| ) could have additional saddle point configurations. However, a quick analysis shows 



that in either case, the full effective potential (periodic in the directions corresponding to 
the gauge field eigenvalues) is attractive both for the scalar and for the gauge field zero 
modes,lii so the eigenvalues are driven towards the clumped saddle point for which the 
KK-mode contributions are negligible. We conclude that the eigenvalues are clumped on 
a scale ( |A.23[ ) whenever 

< — and — — < — . (A.24) 



-R1-R2 / Ri V-R1-R2 / R 

We see that the eigenvalues spread out relative to the sizes of the dual circles as either Ri 
or R2 increases, suggesting a possible phase transition when either of the inequalities in 
( |A.24|) is violated. However, the effective coupling of the lightest KK-modes is Ao/m^^, 



so when the relations ( [A.24| ) are not satisfied, either the Ri or the R2 KK-modes become 
strongly coupled and we can no longer trust our perturbative results for the effective 
potential. Therefore, other methods (discussed in the main text) are required to deduce 
the presence (or not) of a phase transition as the eigenvalues spread. 

74.5. Integrating out massive scalars 

Before closing this section, we describe the result of integrating out a set of p very 
massive adjoint scalar fields in a two dimensional gauge theory on T^. In general, the 
one-loop result is given by ( |5.18| ), where is the covariant derivative for the adjoint 
representation. Here, we consider the special case where the gauge fields are constant 



It is not immediately obvious that the infinite set of image charges does not conspire to 
produce a larger effect, but because of cancellations between rh and —rh terms in the potential, it 
may be checked that the sum of contributions from all charges is of the same order of magnitude 
as the contribution from the nearest charge. 

There is one qualitative difference between the bosonic and supersymmetric cases. In the 
supersymmetric case, the potential between two eigenvalues actually goes to infinity as the eigen- 
values approach antipodal points in the anti-periodic directions, since in this configuration there 
is a (repulsive) image charge from one eigenvalue sitting on top of the other eigenvalue. 
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commuting matrices, so that the holonomies are U = e*^^"^^ and V = e^^^-^^ . Then 

5e//-|hl(dct(-2^J+M2)) 

= |tr(ln(-D? - Dl + M^)) 



\m,n \ ^ ' ^ 

where in the last hne, we have used a result from appendix B. 



Appendix B. Infinite products 

In evaluating the determinants involved in integrating out Kaluza-Klein modes and 
massive scalars in the main text and in appendix A, we encounter infinite products whose 
logarithm gives an infinite sum of the form 

P(5) = ^ln(^2^(m + a)2). (B.l) 

m 

Here, the sum runs over all vectors m in d dimensions with integer components. Now, P 
is clearly periodic in each component of a, with period 1. Thus, we can write 

P{S) = J2 e^'^^^-^Pg, (B.2) 

k 

and we can compute the Fourier transform of P : 



Pk 



-2'Kik-a 



= 1 dai--- [ dad ln(6'^ + (m + af) e 
Jo Jo - 

pmi + l pmd + l 

= J2 dai--- dad H^^ + a^) e'^'^^'^'^ 

^ ''mi Jmd 
= j da ln(^2 ^ ^-2mk-a 

^e-"^' J da e-"^'-^'^*^-" 

„ d f°° da ^a2 TT 

= Coo-7rW ^-je-«^ — 
Jo a^^ 



(B.3) 



= Coo-2^Kd(27r|fc|^). 
Here, Cqo is an infinite constant independent of k and $. 
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Appendix C. Effective Action for the Wilson Line in d = 1 Gauge Theories 

In this appendix, we consider the one dimensional Euclidean gauge theory on a circle 
of circumference R with p scalars and action 

S = j d^tr(^A$^A$^ + Y$^$^-^[$^,$,][$^,$,]). (C.l) 

The analysis here follows that of which considered the special case of two scalar fields. 
Note that here we have rescaled the gauge field and the scalar fields by a factor of g 
compared to our previous analysis, so as to have canonical kinetic terms. We would like 
to integrate out the scalar fields for small values of the dimensionless 't Hooft coupling 
g'^N/M^, to obtain an effective action in terms of the Wilson line of the gauge field around 
the circle. We choose the gauge dtAo = 0, so that Aq is a t-independent Hermitian matrix 
Aq = a. With this choice, the Wilson line is simply given by 

U = e^"^ . (C.2) 

We define 

eM-SeffiU)) = I [ci$,]e-^(*'«), (C.3) 
in terms of which the partition function is given by 

Z = J DUeM-Seff{U)) . (C.4) 

As explained in section 4 of 0, the Haar measure DU arises from the initial measure 
[dAo] upon introducing the Fadeev-Popov determinant associated with the gauge fixing 
condition dtAo = 0. 

By gauge invariance, the effective action must be some function of the variables 

ur. ^ ^tr(C/-) . (C.5) 
At one-loop order, the calculation of S^ff was described in 0, with the result 

oo ^ 

^eVr'(^) = E - P^'')M^ (C.6) 



n 



where x = e . Here, the x-independent term is a rewriting of the Vandermonde 
determinant obtained in writing the Haar measure in terms of eigenvalues. For x < 1/p, the 
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one-loop effective action is positive definite, minimized by tlie saddle-point configuration 
Un = 0. As X passes Xc = the mode ui becomes unstable and condenses to its 
maximum allowed value (as long as all other m„ = 0) ui = 1/2, giving rise to a first order 
large phase transition in the strict g'^N = limit. 

As described in section 6 of , in order to determine the nature of the phase transition 
for weak but non-zero coupling, it is necessary to take into account higher order terms in 
the effective action. The relevant physics may be deduced easily from the effective action 
for wi, obtained by integrating out both the scalar fields and all the modes Wn>i near the 
transition. This takes the general form 



where in perturbation theory b starts at order A^. If the coefficient b is positive at the 
value Xc of x where drops to zero, we will have a second order phase transition with 
the eigenvalue distribution for U changing continuously. On the other hand if (as we will 
find below) the coefficient b is negative at this value x = Xc, the potential develops a 
second minimum which is lower than the first already at some slightly lower value of x, the 
eigenvalue distribution changes discontinuously, and we have a first order phase transition. 
The leading order contribution to b is given by 



where C2, and are the leading coefficients of the terms |w2p5 (^2^-1 + U-2U1), 
and 1^1 1^ in the effective action obtained by integrating out the scalars; these terms first 
arise at one, two, and three-loop order, respectively. We now proceed to compute these 
coefficients, together with the corrections to the coefficient A = N'^m\ of needed to 

determine how the phase transition temperature varies with the coupling constant. 
The higher loop corrections to the effective action are given by 



(C.7) 






2 



2 



(C.9) 



This may be evaluated in perturbation theory using the propagator 



(($z)fcz(il)($j)mn(t2)) = ^knitl -t2,a)Sl 



m 



(C.IO) 
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where the matrix A is defined by 



Aoit 



-tM 



AM 



2M \ 1 — QiaR Y — QM/dgiaR 



(C.ll) 



More details about the perturbative evaluation and our conventions may be found in [15 



Setting M = 1 for now, the contribution to the effective action from the two-loop figure 
eight diagram is given by 



S. 



2 loop _ 9 



/3(p2-p)tr(A(0,aab)A(0,a„c)) 



'eff - 2 

where we have introduced the notation aab = 0(a — c^b^ with 

aa = a®l®l, a6 = l®a(g)l, etc = l®l®a- 
This leads to quadratic terms 

Sl!:7 = —N^Kp' - P) ln(x) J2{x'- + 2xnK\^ 

n 

and leading order cubic terms 

SluZ" = -In'Kp' - P) ln(x)(x2 + 2x'){u2u\ + u.^ul) + ... 



(C.12) 



(C.13) 



(C.14) 



(C.15) 





3a 



3b 



3c 



Figure 21: Three- loop diagrams contributing to the effective action. 

At three loops, the diagrams 3a, 3b, and 3c shown in figure 21 give the following 
contributions 

Slh = -\9\P' -P)l^ j dttr{A{t,aab)A{t,abc)A{t,a^d)A{t,ada)), 



Seff = 1)^/3 / rfttr(A(0,a„b)A(t,Q;„c)A(t,Q!ca)A(0,a„d)), 



^eff 



(C.16) 



(7Xp-1)'/3 / cittr(A(0,a„fa)A(t,a„,)A(t,a,„)A(0,aed)). 
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These contribute to the quadratic term for the lowest mode 

SgLd = -^N''X^\ui\^{p^-p) ln(x)x(2xln(x) - X -3), 

S?uad = -^N^>^^\ui\Mp - 1)' ln(x)x(ln(x)(2x + 1) - 3(x + 1)), (C.17) 



16 
1 

16 

The same diagrams also give the leading |tr([/)|^ terms 



SqLd = -^N^>^^\ui\Mp - 1)' ln(x)a:(ln(x)(x + 1)^ - {x^ +2x + 3)). 



Sfuart = --N'X'\U^MP - l)'ln(x)x3(ln(x) - 1), (C.18) 
SiLrt = -^N^X'\ui\Mp - 1)' In(x)x2(2x2 ln(a;) - x^ - 2). 
Collecting our results, we have 

Seff = A\uif + Bi\u2f + C2(W2W-1 + u-2ul) + Ds\ui\^ + . . . (C.19) 



where 



A/N"^ =(1 -px) - ^A(p2 -p) ln(x)(x^ + 2x) 

A^(p^ — p)xln(x)(2xln(x) + x^ — x — 3) 

16 

- —X'^p(p - l)2xln(x)(ln(x)(x2 + 4x + 2) - x^ - 5x - 6) + C(A^), 
16 

Si/iV2=l(l-px2) + 0(A), (C.20) 
C2/N^ = - lx{p^ - p) ln(x)(x2 + 2x^) + ^(A^), 

o 

D^/N^ = - — A2(p2 - p) lii(x)x2(2 ln(x) + 2x^ - 5) 
32 

- —X^pip - I)2ln(x)x2(ln(x)(2x2 + 3x) - x^ - 3x - 2) + O(X^). 
16 



Prom these above expressions and from (|C.8| ) we find that 

t2 _ 1 (p - 1) ln(p) 
32 p3 



6(x = Xc = l/p)/N^ = -1^— ^^(ln(p)(9p2 + 2p) + V + V - 4p - 4)A2, (C.21) 



which is negative for all p > 1. Thus, the phase transition is of first order at weak coupling. 
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To determine the precise transition point as a function of A, we note that to order A^, 
the second minimum = 1/2 wiU become dominant at the point Xc where 

m?(x)l+6(x)l=0. (C.22) 

Solving perturbatively for Xc as a function of A, we find 



1 , (p-l)(2p + l)ln(p) 

Xc = - + A —: 

p Ap^ 

+ A^ ^)^^(^) ((16p4 _ 33^2 _ ;^Q^x ^ _ gQ^4 _ 20p3 ^ ^^1^2 ^ ^^p + 4) + OiX^). 

128p* 

(C.23) 

Expressing this as an equation for t = l/(i?A^/"^) as a function of m = M/X^^^ gives the 
relation 



Appendix D. Pure and Deformed Yang-Mills Partition Functions on 

In this appendix, we review the exact solution of pure Yang-Mills theory on 



and use techniques developed by Gross and Taylor pl] , p2| , |33[] in order to write the result at 
large in terms of an effective action for one of the two holonomies, U and V. We then 
proceed to study the partition function for a class of deformations which are encountered 
in section 5. 

D.l. Pure Yang-Mills theory on 

An exact expression for the partition function of pure Yang-Mills theory on is 
easily obtained, following Migdal 0], if we use a particular lattice regularization of the 
theory. We put unitary matrices Ul on the links on the lattice, and consider the partition 
function 



^ T D 



L P 

where P denotes the plaquettes of the lattice. Up is a, product of the U^s around the 
plaquette P, and Zp is a plaquette action chosen so that the continuum theory coincides 
with Yang-Mills theory (for instance, the Wilson action Zp = exp(^^tr(t/p + Up^))). 



YM 



Migdal noted that, upon integrating out various links, one is left with an action of the 
form (pj.lD with a plaquette action for the remaining (larger) plaquettes that approaches 



a weighted sum of characters xr^ the representation R of the gauge group : 

Zp{Up) = J2dRXR{Up)e-^^^^^\ (D.2) 



R 
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where the representation R has dimension dji and quadratic Casimir C2(-R), A is the 
dimensionless 't Hooft couphng, and A is the area of the plaquettei^. Using standard 
orthogonahty relations, it is easy to verify the additivity property of (|D.2|) , 

J DUZp{ViU, XAi)Zp>{uW2, XA2) = Zp+p,{ViV2, A(Ai + A2)), (D.3) 

at which point ( p.l| ) can be reduced to an integral over two matrices, U and V , corre- 
sponding to the products of the Uls over the links along the two non-trivial cycles of the 
torus : _ 

Zy^ = I DUDVj2dRe~^^'^''^XR{UVU-^V-^), (D.4) 

R 

where A is the 't Hooft coupling times the area of the T^. We can then use the result 

J DUxR{UAU-'B) = j-XR{A)xR{B) (D.5) 
to integrate out one matrix, say U, from ( p.4| ), obtaining 

Zym = [ DV $^e-^^^(^)xfi«i^(V^)- (D.6) 

R 

This integral is also trivial to perform, since R ^ R contains the identity precisely once, 
yielding the well-known expression for Zym purely as a function of A, 

^.rn = J]e-^^^(^). (D.7) 

R 

In the rest of this subsection, we shall be interested in determining, in the limit of large 
A^, a more explicit form for the effective action, Seff{V), corresponding to the partition 
function ( |D.6D . 



Thus, we seek a large expansion of the quadratic Casimir and character appearing 
there. As pointed out by Gross and Taylor in pT| - p3| , one may obtain the correct large 
N expansion of (p.6| ) by summing independently over representations formed from tensor 



20 



In Q, Migdal considered a more general situation in which ( p.l ) was not required to yield 



pure Yang-Mihs theory in the continuum hmit. As such, he originahy wrote the factor in the 
exponent in terms of a representation-dependent couphng, qr, noting that gj^ gy j^jC2{R) gives 
the correct Coulomb law. One can also see this directly by applying Migdal's analysis to the 
Wilson action. 
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products of finitely many fundamentals and anti-fundamentals. The Young tableaux for 
any such representation may be obtained uniquely by adjoining the Young tableau for 
some representation obtained from the product of finitely many anti-fundamentals with 
the Young tableau for some representation obtained from the product of finitely many 



fundamentals. Following the notation in [|31| , |33|] we denote these representations by S and 
R respectively, and denote the representation corresponding to the original tableau, termed 
the "composite" representation of R and S", by T = SR. 

Letting ur (ns) denote the number of boxes associated to the representation R (S), 
we may express the quadratic Casimir of the composite representation T in terms of those 
of its components as: 

C2(T) = C2{R) + C^iS) + (D.8) 

where 

C2{R)=NnR + 0{N'>) , C2iS) = Nns + OiN"^) . (D.9) 

Denoting the set of Young tableaux with n boxes by T„, and denoting both an arbitrary 
element of T„ and the representation corresponding to it by & T^, we may now write 
( PXel ) at large as 

n,n'Y„eT^Y',eT^, 

n' 

We now seek to expand the characters appearing in the integrand as a polynomial in 
traces. To proceed, we first define some further notation. To any given tableau, G T„, 
we may associate not only an irreducible representation, also denoted by Y^, of SU{N), 
but also an irreducible representation, Y-a, of the permutation group Sn- We use Xy„{V) 
to denote the character of the SU{N) representation Y^ evaluated on V E SU{N), and 
Xr; (cr) to denote the character of the 5"^ representation Y^ evaluated on a G 5"^. To each 
tableau, G T^, we also associate a conjugacy class, of Sn, which can also be labeled 
by the numbers {cTj}, which specify the number of cycles of length i in an element of py„. 
Finally, we let T^r denote the Schur function associated to a G 5",^, 

rAV)= XyJ^)XYjV)=l[tr{Vr'. (D.ll) 

The character of any representation Y^ G is easily expanded in terms of Schur functions 

^>'"(^) = ^ E XyJ^^AV). (D.12) 

o-eSn 
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To study ( p.6|) , we shall need the following generalized version of (p.l2| ) for composite 



representations that was worked out by Gross and Taylor |33 



XK,y^^) = ^ E E XyJ-)Xy,Srn.AV.vh (D.13) 



where 



I 

min(cr;,r;) / \ / \ 

T,(o.(o {V, yt) = J2 ^j^) [I) {-l)H'kltv{vY^-hv{V- 



(D.14) 



l\Ti-k 



Using this result, we may now write ( |D.6| ) as 



«-("+»■)* X 



z,„.= / DvY. E E E E 

n,n' <t,<t'6S„ t,t'6S„/ y„eT„ y^,eT„, 

X9 K)X9' (^') ■ . 



(n!)2(n'!)2 

(D.15) 

The sum over Young tableaux can be performed using the completeness relation 

where p^r denotes the conjugacy class to which a belongs and |pcr| its dimension. Using 
this result we obtain 



^^EE E T,.(y,yt)T,,(y,yt) 

/ ^ o ^ e lt.lt . 

n,n' crEbn TEh^i 



X 

n,\n!\ J-J.'-L ^ ' ^ 'J 

n,n' aESn tES^ 

min((T;,r;) 



(D.17) 



A;,A;'=0 
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We now replace the sums over cr, r with (appropriately weighted) sums over the ai.Tj 
which label conjugacy classes. Typically, such sums must be subject to the restrictions 
Ylii "^^i = ^5 Ylij j^j — ^' are difficult to evaluate. Fortunately, the sums over n, n' 
serve to lift these restrictions, leaving us with the following expression 



oo min((T;,r;) / _ Ai. 



X ■ 



I \ o-i,Ti=0 k,k'=0 



I 



k\k'\{ai - k)\{ai - k')\{Ti - k)\{Ti - k')\ 
which can be rewritten as 



[tT{V')t^{V-')Y'^''' 
' [tr(yOtr(y-^)]"^'+''^ 



(D.18) 



Zym= DV\[l E 



<T+T+2 max(A;,fc') 



I a,T=0 k,k'=0 

a + max(/c, k')]l [r + max(/c, k')]l 
[cr+ |A;-A;'|]![r+ |A;- A;'|]! 



aWlklk'l 



[tr(y')tr(y-')] 



(D.19) 



We now consider the sum at fixed /, which may be rewritten, after a little algebra, as 



oo oo 
cr,T,fc'=0 fc= — oo 



1 



e 2 



cr + T 



-XI 



aWlk'l I 



-XI 



X [tr(yOtr(y-0]'^+"+l'l 
In our large N limit, ( p.2Cl| ) becomes 



{a + \k\ + k'y.jr + \k\ + k'y: 

(a+|A;|)!(r+|A;|)!(|A;| + A;')!. 

(D.20) 



exp 



2e~^ - e~^^ 
I 



tr(y^)tr(y-') 



(D.21) 



with corrections that are subleading in in the exponent. As a result we find that ( p.l7| ) 
can be written as 

Zy^= [ DV exp\-Seff{V,X)}, (D.22) 



where 



Seff{VA) = Yl 



1=1 I 



e~- - 1 - 1 



tr{V%{V- 
I 



(D.23) 
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To study this result, we write it in terms of the moments = tr{V"')/N of the 
eigenvalue distribution, 



. 2 

— 1 



\Ur 



|2 



n 

n=l 



> . (D.24) 



A model similar to this was studied in 0. Here, the masses of all moments, Un, are positive 
for all positive values of the coupling. A, and thus the dominant eigenvalue distribution 
is the uniform one with Un = for all nonzero n. At large A^, the Un are essentially 
independent variables and can be integrated out to yield 



-2 

A n 



Zyn.--llil-e-^j . (D.25) 

This is in agreement with the known partition function for pure Yang-Mills theory on 
at large N. 

D.2. Deformed Yang-Mills theory on 

In this appendix, we derive a formula relevant to our study in section 5 of pure Yang- 
Mills theory deformed by very massive scalars. In section 5.5, we argue that in a certain 
regime of parameter space, integrating out very massive adjoint scalars gives, to a good 
approximation, an expression ( ^.24|) for the partition function, where / and g depend only 
on the eigenvalues of U and V respectively. Now, changing variables U WUW~^ and 
integrating over W (which doesn't change the result since the integrand cannot depend on 
W) we obtain 

Zym = f DW DU DV YdRe-^^^^^^XR{WUW-^VWU-^W-^V-^)e-f^^^e-3^^\ 

(D.26) 

Here, we have used the invariance of both the measure DU and of the function f{U) under 
the transformation used in the change of variables. To proceed further, we seek to evaluate 
the integral 

I{Di,D2) = j DCxr{CDiC-^D2CD^^C-'D^^). (D.27) 
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Using SU{N) symmetry, the form of /(-Di, D2) is restricted to 



(D.28) 

To compute /(Di, 1^2)5 we thus need only to determine the three numbers a, /3, 7. To 
constrain their values, let us first look at /(I, 02). Equating ( p.27|) and (|D.28|) we obtain 

a + P[dl + xr{D2)xr{D^')] + ldlxR{D2)XR{D2^) = dR. (D.29) 

We next consider integrating J DDiI{Di, D2). Comparing the result obtained by inter- 
changing the integrations over Di and C yields 

a + P[l + xr{D2)xr{D2')] + -iXR{D2)xR{D2^) = ^Xr{D2)xr{D2^) . (D.30) 
Equations ( p.29| ) and ( p.30| ) give 4 equations relating the coefficients a, /3, 7 : 

a + (3d\> = dR, 



'R 

/3 + 74 = 0, 
a + p = 0, 

dR 



This system has a unique solution 



d 



a 



R 



dl-1 



4-1 



1 



7 



(D.31) 



dR{dl-iy 

Applying ( p.28|) with these values to our expression ( p.26| ), we may rewrite Zym as ( |5.25|) . 



Appendix E. More about the Monte-Carlo simulations 

The Monte-Carlo simulations presented here were written using an elementary im- 
plementation of the Metropolis algorithm. The ensembles were fully thermalized between 
samplings with various time-time correlators being checked. 
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In the case of the 0+0 matrix integrals the implementation is extremely simple. For the 
0+1 matrix quantum mechanics we require a lattice of L spatial sites, with circle topology, 
and each site is equipped with the scalar adjoint matter matrices, with the gauge field 
living as usual on the links. Since the gauge dynamics is trivial in this low dimension 
we may perform a gauge transformation to make the unitary link variables equal on all 
sites. Furthermore we may use up the remaining gauge freedom by diagonalizing this 
unitary link matrix. The remaining diagonal components are pure phases, and are physical, 
giving the eigenvalues of the Polyakov loop when raised to the power L, the number of 
lattice sites. Wc carefully ensure that the Jacobian introduced by this gauge fixing is 
properly implemented. This unitary matrix measure is the non-perturbative version of the 
Vandermonde determinant and is given by sin^ ^{^i ^ ^j) for i = 1, ■ ■ ■ , A^, where 
^lOi g^j,g ^Yie eigenvalues of the unitary link. This measure factor is implemented by taking 
its log and introducing it as a potential term in the action. We automatically adjust the 
Metropolis step size to ensure decent acceptance rates, and use independent step sizes 
for the unitary link eigenvalues, and for both the scalar adjoint matter diagonal, and 
ofi-diagonal components. 

For the large N behaviour we study we require relatively few lattice points to accu- 
rately capture the continuum behaviour. For p scalars with p = 2,4 the data presented 
here uses 10 lattice points. For the p = 9 data 5 lattice sites were used. For individual 
values of A and M we checked that this was sufficient for the quantities we measured, 
finding that doubling or quadrupling the number of lattice sites did not change the results 
at the level of one percent. 
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